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Future or ongoing galaxy redshift surveys can put stringent constraints on neutrinos masses via 
the high- precision measurements of galaxy power spectrum, when combined with cosmic microwave 
background (CMB) information. In this paper we develop a method to model galaxy power spectrum 
in the weakly nonlinear regime for a mixed dark matter (CDM plus finite-mass neutrinos) model, 
based on perturbation theory (PT) whose validity is well tested by simulations for a CDM model. 
In doing this we carefully study various aspects of the nonlinear clustering (nonlinear neutrino 
perturbations and the higher-order growth functions), and then arrive at a useful approximation 
allowing for a quick computation of the nonlinear power spectrum as in the CDM case. The nonlinear 
galaxy bias is also included in a self-consistent manner within the PT framework. Thus the use of 
our PT model can give a more robust understanding of the measured galaxy power spectrum as 
well as allow for higher sensitivity to neutrino masses due to the gain of Fourier modes beyond the 
linear regime. Based on the Fisher matrix formalism, we find that BOSS or Stage-Ill type survey, 
when combined with Planck CMB information, gives a precision of total neutrino mass constraint, 
er(m„,tot) — 0.1 eV, while Stage-IV type survey may achieve er(m„ lto t) — 0.05 eV, i.e. more than a 1- 
a detection of neutrino masses. We also discuss possible systematic errors on dark energy parameters 
caused by the neutrino mass uncertainty. The significant correlation between neutrino mass and 
dark energy parameters is found, if the information on power spectrum amplitude is included. More 
importantly, for Stage-IV type survey, a best-fit dark energy model may be biased and falsely away 
from the underlying true model by more than the 1-cr statistical errors, if neutrino mass is ignored 
in the model fitting. 

PACS numbers: 98.80.Es,14.60.Pq,98.65.Dx 



I. INTRODUCTION 



The concordance A-dominated cold dark matter (ACDM) model for structure formation in the Universe is remark- 
ably successful in describing various data sets such as cosmic microwave background (CMB) anisotropies, Type-la 
supernova distance measurements, observations of galaxy clustering and cluster counts, and weak gravitational lens- 
ing (e.g., [l|; 0)- However, the concordance model requires that the present-day energy budget of the Universe is 
dominated by unknown two dark components. One is dark matter that is needed to explain the nonlinear aspects of 
gravitational clustering in structure formation, and the other is the cosmological constant contribution or perhaps a 
more generalized form dubbed as dark energy, which is needed to explain the cosmic accelerating expansion. Under- 
standing the nature of these dark components is one of the most important, profound problems in cosmology as well 
as particle physics. 

We now know that the Big-Bang relic neutrinos contribute to dark matter energy density by some small fraction, 
because the neutrino oscillation experiments @, 0, H, Q have shown that neutrinos have finite masses (also see @, H| 
for a thorough review). However, the oscillation experiments are sensitive only to mass square differences between 
different flavor neutrinos, therefore the most fundamental constant of neutrinos, absolute mass scale, is not yet known, 
although the lower bound on total neutrino mass can be inferred as m„.tot ~ 0.06 or 0.1 eV for the normal and inverted 
mass hierarchies, respectively. On the other hand, the direct experiment has put only a weak upper bound on electron 
neutrino mass such as m„ c ^ 2 eV . 

Cosmological probes can give a complementary, albeit indirect, method for constraining neutrino masses. There are 
two kinds of the methods. First is via the effect on cosmic expansion history. If neutrino species are massive enough 
as m v ^ 0.5 eV, the neutrinos became non-relativistic before recombination epoch and then imprint characteristic 
signatures onto structures of the CMB anisotropy spectra [l(3, El- On the other hand, low-redshift geometrical 
probes such as Type-la supernovae (e.g. [12]) and the Baryon Acoustic Oscillation (BAO) [l3| are sensitive to the 
present-day energy density of non-relativistic matter (f2 m o) that is given by the sum of CDM, baryon and neutrino 
contributions: f2 m o = f^cdmo + ^bo + ^vo- Therefore combining these geometrical measurements can constrain 
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neutrino mass: for example, [l[ already succeeded in obtaining the presumably best-available constraint from this 

method, m„,tot ~ 0.6 cV (95% C.L.). 

Alternative cosmological method is using clustering information of large-scale structure. Due to large ve- 
locities of the frozen Fermi-Dirac distribution, neutrinos cannot cluster on scales below the neutrino free- 
streaming scale that has a characteristic dependence on neutrino mass and redshift as given by fcf Sj i ~ 
0.023/iMpc _1 (m Z y i i/0.1eV)(r2 m o/0.23)[2/(l + z)} 1 / 2 , comparable with the BAO scales for neutrino mass scales of 
interest. As a result, the presence of finite-mass neutrinos suppresses the amplitude of low-redshift power spectrum 
on the small scales by at least the amount of a few percent, compared to the model without finite-mass neutrinos, 
for a fixed f2 m o [13]. Thus given the CMB normalization of primordial power spectrum, total neutrino mass can be 
explored by measuring clustering strengths of low-redshift large-scale structure via galaxy redshift survey [T5|.[l6l [ItJ , 
weak gravitational lensing [l8l. Lyman-a forest power spectrum pol ] and potentially 21cm observations j2l|7 The 
existing data sets have put more stringent upper bounds on neutrino mass, m u ,tot ~ 0.2-0.5 eV, than the direct 
experiment limit, although some residual systcmatics are under discussion. 

There are a number of ongoing and planned galaxy redshift surveys such as WiggleZ [97|, FMOS [2^], BOSS l9g| . 
Subaru redshift survey known as the former project WFMOS [H, HETDEX [To^, EUCLID [ToJ, and JDEM [l02j j. 
The primary scientific target of these surveys is exploring the nature of dark energy via the BAO experiment. At the 
same time these surveys promise to achieve the high-precision measurements of galaxy power spectrum amplitudes to a 
percent level precision at each wavenumber bins, and therefore offer a possibility to dramatically improve cosmological 
constraints including neutrino masses [23|, [24|, [25[ . 

Thus large-scale structure probes are very promising, however, the main obstacle is nonlinear effects such as 
nonlinear gravitational clustering, galaxy bias and redshift distortion. Recent theoretical studies have shown 
that, even at scales as large as ~ 150/i _1 Mpc relevant for both the BAO and neutrino free-streaming scale, 
the standard linear theory, which gives remarkably successful agreement with CMB measurements, ceases to be 
accurate. The nonlinear effects are found to be sig nificant compared to the precision of future surveys, us- 
ing A-body simulations p6|. l27l. Epj. [29l. l30l . l3ll . l32j and analytical studies inspired from perturbation theory 
[33l . l34l . [35l . [H, [37], [H, [391 401411. |42i |4SI|. However in most of these studies the contribution of finite-mass neu- 
trinos were ignored. The nonlinear effect of finite-mass neutrinos on the power spectrum needs to be understood in 
order to attain the full potential of future surveys, which is also important to minimize the possible systematic error 
on BAO experiments caused by the incorrect assumption that neutrinos are massless. 

Therefore the aim of this paper is developing a formulation to model nonlinear galaxy power spectrum in a mixed 
dark matter (CDM plus finite-mass neutrinos) model, based on standard perturbation theory (SPT) (see [45[ for 
a thorough review of perturbation theory for a CDM model). Here we mean by "standard" that the next-order 
corrections to the power spectrum, i.e. the one-loop corrections, are included. In doing this we carefully study 
various aspects of the nonlinear clustering: estimate the nonlinear neutrino perturbations by solving the collision- 
less Boltzmann equation hierarchies and study the higher-order growth functions of CDM plus baryon perturbations 
that have complicated scale- and redshift-dcpcndcnccs similarly to the linear-order growth rate. Then, given the 
detailed assessment of various effects, we will arrive at a useful approximation to compute the nonlinear matter power 
spectrum whose results were highlighted in [|(| . We then include a modeling of nonlinear galaxy bias self-consistently 
within the SPT framework following the method developed in (47[. Thus, while the nonlinear redshift distortion effect 
is not yet included, our model of the galaxy power spectrum can be compared to the actual measurement such as 
that in [T3], where the redshift distortion effect is removed using the Finger-of-God compression algorithm [48|. For 
preparation of such a study we will demonstrate parameter forecasts for neutrino mass constraints expected from 
some of galaxy surveys listed above, paying a particular attention on the correlation between neutrino mass and dark 
energy parameters in the weakly nonlinear regime. We also discuss a possible systematic error in the constraints on 
the dark energy parameter caused by the neutrino mass uncertainty. 

The structure of this paper is as follows. In § [IT] we develop the formulation of SPT method for computing the 
nonlinear matter power spectrum. In § Mil we then study the effect of finite-mass neutrinos on the matter power 
spectrum by varying the neutrino masses within the range inferred from the constraints. After including a model 
of nonlinear galaxy bias based on perturbation theory in § IIV1 we study parameter forecasts of neutrino masses and 
dark energy parameters using the Fisher matrix formalism in § [Vj § I VII is devoted to summary and discussion. Unless 
explicitly stated, throughout this paper we assume the concordance ACDM-like cosmology with finite-mass neutrino 
contribution, which is consistent with the WMAP results [l|. The fiducial model is: the density parameters are 
f! m Q = 0.24, £l m oh 2 = 0.1277, and fibo^ 2 = 0.0223. The neutrino effect is studied by varying the neutrino mass 
scale. For simplicity the number of neutrino species is assumed to be N v = 3 because the matter power spectrum is 
sensitive to the sum of neutrino masses, m„ it ot = N v m v . We assume a flat universe and consider wq = — 1 for dark 
energy equation of state. For the primordial fluctuation parameters, the amplitude, the tilt, and the running, are set 
to = 2.35 x 10 -9 , ns — 1.0, and as = 0, respectively. 
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II. PERTURBATION THEORY FOR NONLINEAR MATTER POWER SPECTRUM IN A MDM 

MODEL 

A. Preliminaries 

First we write down basic equations to describe structure formation in a MDM model. Throughout this paper, we 
focus on the evolution of matter fluctuations consisting of MDM (CDM plus massive neutrinos) and baryon: 

_ 5p c + Sp h + 5p v 

o m = = JciAb + Jvd v , (1) 

Pm 

where the subscript 'm', 'c', 'b', V and 'cb' stand for total matter, CDM, baryon, massive neutrinos, and CDM plus 
baryon, respectively, and <5 c b and 5 V denote their density perturbations. The coefficients, / c b and f v , are the fractional 
contributions of each component to the present-day total matter density: 

r _ ^cO + ^bO p _ ^vO _ 1 t r^i m v,tot ,cy\ 

Jcb — p: > Ju — — 1 — /cb — TTTTri U2 > \ ' 

with the density parameter, f2i0j being defined as fi^o = SttG Pi(to) / (3Hq) (i = m, c,b, u) where the parameter, h, is 
dimcnsionlcss Hubble constant defined as Hq = lOOh kms _1 Mpc~ 1 . In the limit of f c h — ► 1, the results shown below 
recover a CDM model which does not contain massive neutrinos. The evolution of homogeneous and isotropic universe 
is controlled by CDM, baryon, massive neutrinos and dark energy whose equation of state is simply assumed to be 
constant in time: £>de = woPde where wq is referred to as the equation of state parameter. Then, the background 
Friedman equations become 

H 2 = — (Pm+PDE), (3) 

H = -'±H\l + w ), (4) 

where the dot ' denotes the derivative with respect to cosmic time, t: a = da/dt, and the Hubble parameter, H, is 
defined as H = a/ a. 

We are specifically concerned with the nonlinear growth of matter perturbations, <5 m . Let us first consider the 
contribution of CDM plus baryon perturbations, S c h, to the total matter perturbations. In order to evaluate the 
nonlinear evolution of <5 c b, we treat the CDM plus baryon components as a single pressure-less fluid. The continuity 
equation and the Euler equation for the CDM plus baryon fluctuations in Fourier space are given in [45|, [4^] as 

g- ^MM) + 9 ch (k;t) = - J 0^ a (k',k-k')6 ch (k-k';t)e ch (k';t), (5) 



1 f d 3 k' 

2 / (2^)3 /?(fc '' fc ~ k ' ] ° ch[k ~ fc ' ; t] 0ch{k ' ; *>' (6) 



where the velocity divergence # c b is defined as # c b = V • v c b/(aH) in real space. Note that we assume an irrotational 
flow, i.e. the vorticity is neglected (50j . The Fourier kernels to describe the nonlinear mode coupling, a and (3, are 
defined as 

fu u-\-t i fcl ' fc2 atu i- \ - (fci -fc2)|fci+fc 2 | 2 
a(feljfe2 ) = ! + __, p( ku k 2 ) = . (7) 

Taking the time-derivative of Eq. §5§ and also using the Euler equation © yield the second-order differential equation 
for 5 ch : 

Sl h + 2HS' ch - '^H 2 {1 - n DE ) 5 m 

= - J ^3 [ a(k', k - fc'){ [H 6 ch (k - k')9 ch (k')\ ■ + 2H 2 S ch (k - k')9 ch (k') } 

+ ^H 2 /3(k',k-k / )9 ch (k-k / )9 ch (k') . (8) 
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Thus Eq. {8j contains S m = / c b<5cb + fv$v and cannot be solved unless the neutrino fluctuation field, <5„, is specified. 

So let us move on to discussion on the neutrino perturbations. Unlike CDM and baryon, the finite-mass neutri- 
nos have a large velocity dispersion following the frozen Fermi-Dirac distribution, and cannot be treated as fluids. 
Therefore, exactly speaking, it is necessary to solve Eq. coupled with the collision-less Boltzmann equations for 
neutrino perturbations that include the nonlinear terms. This is still computationally expense, especially for solving 
the nonlinear Boltzmann equations. Here we rather consider the approximated method for solving the nonlinear 
perturbations as described in the next subsection, and will also assess an accuracy of the approximation. 

In our method we focus on the linearized collision-less Boltzmann equations for neutrino perturbations [5l| : 

% = -^ 1+H /-±I°, (9) 
at a In q 

, qk , . tk , d In fn , . 

*i = :M*o - 2*2) - — cp— A (10) 



Sat 3aq dlnq' 

Qk ■ 
(2£+ l)ae l 



% = , n , q \ n [^-i -V+ W+i] (^>2), (11) 



where the variables, q and e, are comoving 3-momentum and proper energy defined as e = (q 2 +a 2 m 2 i ) 1 ^ 2 , respectively, 
and the function <j> is the gravitational potential perturbation under the conformal Newtonian gauge (see below). The 
superscript ' denotes the derivative with respect to conformal time. The function fo is the zeroth-order (isotropic) 
Fermi-Dirac distribution, given as /o = 2/(e £ / aT + 1), and the function 'J is the linear-order perturbed distribution. 
The full phase-space distribution function of neutrinos is given in the linear regime as 



f(x\ qj /a, t) = / (g)[l + *(x*,q, hj,t)], (12) 

where the momentum vector is rewritten as qj — qfij with hjh J = 1. The variables, 4'^, appearing in the Boltzmann 
equations above are the ^-th moments in the Legendre expansion of ^: 

oo 

*(fc, n, q, t) - ^(-») f (21 + l)9 t (k, q, t)P e (k ■ n), (13) 

e=o 

where Pi is the Z-th order Legendre polynomial. The neutrino density perturbation is given by integrating the 
monopolc contribution of neutrino perturbations over momentum: 

MM) = f q 2 dqef Q (q)^ Q (k iq ,t). (14) 

The system of momentum hierarchies, Eqs. ©, ([TU| . and (fTTj) . can be solved once the gravitational potential cf> 
is given. One of the Einstein equations, the Poisson equation, relates the potential <f> to the total matter density 
perturbation S m on subhorizon scales: 

- k 2 <f){k, t) = AirGa 2 p m S m - (15) 

On scales smaller than the neutrinos' free-streaming scale, k ^ fcf s , the neutrino perturbation would be absent, and 
the Poisson equation roughly becomes —k 2 4>(k,t) w 4TrGa 2 p m f c bS c b- Thus on these small scales the dynamics of 
neutrino perturbations are governed by the CDM plus baryon perturbations. We have so far written down all the 
basic equations that govern the dynamics of density perturbations for each components, <5 c b and b v . A quantity that 
is more relevant for actual large-scale structure probes such as galaxy clustering is the power spectrum of total matter 
including nonlinear corrections: 

(Mfe;t)Mfc';i)> = (27r) 3 S D (k + k')P ni (k;t). (16) 

The power spectrum, P m , is defined in terms of the density perturbations of CDM, baryon and neutrino perturbations 

as 

P m (k; t) - f 2 h P ch (k; t) + 2/ cb /„ P chu (k; t) + f 2 P v {k; t), (17) 
where P c bv{k) is the cross spectrum between 5 c b and <5„. 
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B. On the treatment of neutrino perturbation 



Strictly speaking, in order to compute the total matter perturbation, S m , in the nonlinear regime, we need to solve 
Eq. ijSJ) coupled with nonlinear collision-less Boltzmann equations for massive neutrinos, which seems computationally 
expensive. In order to avoid this obstacle, in this paper we employ a simple approximation that allows to analytically 
compute the nonlinear power spectrum in a MDM model based on the standard perturbation theory (SPT) (see [46| 
for the similar discussion). 

Let us begin with recalling characteristic properties of neutrino clustering on scales up to k ~ 0.1 ftMpc - . Firstly, 
the neutrino perturbations contribute to nonlinear gravitational clustering via its contribution to the gravitational 
potential, where, implied in Eq. (|15p . the perturbation of physical neutrino density, 5p u = p u 5 u , affects the gravitational 
potential. Thus the contribution is suppressed by a small factor f v , currently limited as f v ^$ 0.05 [l[ , even if the density 
perturbations of CDM and neutrinos are in similar amplitudes as predicted by the adiabatic structure formation 
scenario. Secondly, the neutrino perturbations would tend to stay in the linear regime due to the large velocity 
dispersion, at least on scales relevant for the BAO scales. These facts suggest that the nonlinear power spectrum 
arises mainly from the nonlinear perturbations of CDM plus baryon. Thus we model the nonlinear matter spectrum 
based on SPT (see below), but including only the linear-order perturbations of neutrinos: 

P* L (fc; t) = P c N b L (fc; t) + 2/ cb /„ i* b „(fc; t) + ft P^(k; t), (18) 

where the spectra with superscript "NL" denote the nonlinear spectra described below, and the spectra with "L" 
are the linear-order spectra. With this assumption, the neutrino perturbations can be precisely computed for given 
initial conditions by using the publicly available codes such as CMBFAST [52[ and CAMB [53[. The validity of our 
assumption is studied in detail in Appendix [SJ Here we briefly summarize the result. As explained around Eqs. (|lip 
and (|15[) . nonlinear clustering of neutrino perturbations is driven by nonlinear gravitational potential supported 
by CDM plus baryon perturbations, in a CDM dominated structure formation model. Therefore, the nonlinear 
correction to neutrino perturbations can be qualitatively estimated by solving the linearized Boltzmann equations 
(|5|)- ((TT|) . where the nonlinear gravitational potential due to the total matter density perturbations given by Eq. (TPS)) 
is inserted into the gravitational force term (ignoring the nonlinear neutrino perturbations). The results are shown in 
Fig. 1101 Nonlinear clustering indeed causes a nonlinear evolution of neutrino perturbations, deviating from the linear 
theory prediction. The nonlinear effect causes greater amplitudes of the neutrino perturbations on larger k and at 
lower redshifts; e.g. the fractional difference between the linear and non-linear density perturbations 8^/8^ reaches 
to ~ 10% on k ^ 1 /iMpc -1 at z — for f v = 0.05. However, the nonlinear effect on the total matter power spectrum 
is suppressed by additional small factor /„ as implied in Eq. (fT8|) . In conclusion the nonlinear correction to the total 
matter power spectrum is smaller than one percent level in the amplitude for a range of neutrino masses, /„ ^ 0.05. 
For these reasons, throughout this paper, we employ the approximation (|18p . where the neutrinos affect nonlinear 
power spectrum of total matter via the effect on the growth rates of CDM plus baryon perturbations as described in 
the next section. 



C. Perturbation Theory Approach 

In this subsection we develop a method to compute nonlinear power spectrum of CDM plus baryon perturbations, 
P c b(k;t) in a MDM model based on perturbation theory. First, in order to solve Eq. J5]), we expand the density and 
velocity perturbations in a perturbativc manner: 

* cb = *$? + « + 4? + • • ■ . = »2? + 0$ + ^ + ■ ■ ■ ■ (19) 

where the superscript '(i)' denotes the i-th order perturbation. Here, we include the next-to-leading order corrections 
for P c b (£:;£), which are expressed as 

P cb (fc; t) = P c L b (fc; t) + P c ( b 13) (k; t) + P c ( b 22) (k; t). (20) 

The first term P^ h denotes the linear power spectrum of CDM plus baryon. The last two terms describe the nonlinear 
corrections, the so-called one-loop corrections, and the superscript '(13)' and '(22)' denote the multiplied order of 

(l) (3) f2) (2.) 

perturbations, (^b "cb ) an< ^ \^cb "cb )• ^ c ^ nus mcm dc contributions up to the third-order perturbations. 

Inserting the formal solutions p^|) into Eq. © gives, at the lowest order of perturbations, the differential equation 

S: h {1) + 2HS' ch {1) - ^H 2 (l - n DE ) <5« = 0. (21) 
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This equation can be straightforwardly solved, together with the linearized Boltzmann equation for neutrino pertur- 
bations (e.g., [IH). In this paper we use CAMB [53[ to obtain the accurate solution of 5^. Before moving on to the 
higher-order perturbations of <5 c b, for convenience of our discussion, we formally write down the linear-order solutions 
of density and velocity perturbations expressed as 0, EH : 



D cb (M)A(ft), 0^(k;t) 



dD ch (k;t) 
din a 



A(fc), 



(22) 



where the quantity A(fc) represents the initial perturbation variables at an early epoch ii n ;, sufficiently in the linear 
regime, e.g. the Compton-drag epoch. The ensemble average gives the initial power spectrum: 



(A(fc)A(fc')) = (2ir) 3 6 D (k + k') ^cb(Mini)- 



(23) 



The effect of massive neutrinos can thus be described as the scale-dependent growth function, _D c b(fc; t). At wavenum- 
bers smaller than the neutrino free-streaming scale, kf s , the neutrinos can cluster together with CDM and baryon. On 
the other hand, at k > fcf s , the growth of CDM plus baryon perturbations is suppressed due to the weaker gravitational 
force caused by the lack of neutrino perturbations. Thus the growth function has asymptotic behaviors given in pij 
as 



D ch (k;t) °c 



Dx(t) (k<k la ) 
D^t) 1 ^ (fc»fc fs ) 



(24) 



where D\(z) is the growth rate for a CDM model without massive neutrinos, but with the same matter density f2 m o 
to that of the MDM model (in this case the growth rate has no scale dependence), and the parameter p is defined as 
p = (5 - V25 - 24/„)/4. 

We now consider the second- and third-order perturbations. Substituting the linear solutions into the r.h.s of Eq. ((SJ) 

(2) 

yields the differential equation for the second-order perturbation S\ : 



x" (2) i oh-a' (2) 
<Jcb + 2Hd ch 



^De) fcbS, 



(2) 



d 3 k 1 d s k 2 
(2tt) 3 



5 D {k-k 1 -k 2 )A(k 1 )A(k 2 ) 



i(fei,fc 2 )| 



d\na 



OIJ 2 dDcb(h) n 
An Ti ^cb 

din a 



(fc)} 



■P{ki,k 2 ) -H — — 

2 a In a a In a 



(25) 



where, as described in the preceding section, we have ignored the second-order contribution of neutrino perturbations, 

/cb<5ct? m deriving the equation above. The formal solutions of Eq. (|25|) 



i.e. set 5^ = 0, and therefore used 6m 
can be written as 



d 3 kid 3 k 2 
(2tt) 3 



^(fe-fei -fe 2 )A(fei)A(fe 2 



a{k 1 ,k 2 )Af ) (k 1 ,k 2 ;t) + !3{k 1 M)Bf\k 1 ,k 2 - 1) 



(26) 



where and Bg are the second-order growth functions given in detail in Appendix IB1 

There arc notable differences between the second-order perturbations in models with and without massive neutrinos. 

(2) (2) 

Firstly, the second-order growth functions A s and B s are scale-dependent originating from the scale-dependence 
of the lineai' growth rate. Thus additional nonlinear mode coupling arises via the scale-dependent growth rate, in 
addition to via the shape of the input linear power spectrum. Secondly, the gravitational force is weaker in a MDM 

(2) (2) 

model because we ignored the second-order neutrino perturbations in our method, i.e., S m = fcbS ch with / c b < 1. 
These imply that the second-order density perturbations are suppressed compare to those of CDM model with same 
^mo- We will in detail show the results below. 

Similarly, a formal solution of the third-order perturbation can be expressed as 



d^k^k^kz 
(2tt)" 
r a ( 3 ) 



S D (k - Jbi - fca - fc 3 )A(fci)A(fe 2 )A(fc 3 ) 

(t)} - a23 I i{a a , 3 cg 2i3 (*) +/?2,3£>fi,3(*)} 

(3) 



2,3$ 
(3) 



-Pi,23{a 2 , 3 E sl23 (t) + f3 2 , 3 F sl , j 



(*)}, 



(27) 
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where the Fourier kernels a(ki, k 2 ), /3(k±, k 2 ) and a(k\, k 2 + £3) are abbreviated as 0:1,2, 01,2 and ai.23 respectively, 

and the third-order growth functions I s (X = A,B,C,D,E,F), abbreviated as X\ (ki,k 2 ,ka) = X$ 123 and so on, 
are given in Appendix [Bj 



D. One- loop Corrections to P c b 



We now study the higher-order density perturbations for a given MDM model. Using the formal solutions, we can 



derive the explicit expressions for the one- loop corrections to P cb . First let us consider P cb '. Using Eq. (f2"o) . the 

(2) (2) / (2,2,) 

ensemble average (<^ b (k;t)5 ch ' (k ;tj) yields the power spectrum P cb . 



P C (22) (M) = 2 J l^ar^y*) +0 1 , 2 B^(t)] 2 P^(k 1 AnOP^(k2AnO 



k — k\ -\-k/2 



(28) 



(2') (2) 

where we have used the abbreviated expressions such as A y s (fc 1; k 2 ,t) = A y s [ 2 (t), and the /^-integration has to be 
done under the condition k = k\ + k 2 . The prefactor 2 arises from the Wick's theorem in evaluating the ensemble 
average, (A(fci)A(fe2)A(fe3)A(fc4)). Changing the integration variables to r = k\jk and /i = k ■ kijikkx), the 

expression of Pj 22 "* is rewritten as 



jL3 poo /*1 

^b 2 MDM(M) = 2^j o r 2 drPt h (kr;t) J d^ P&(ky/l + r 2 - 2r/i; t)JCf (k, r, m t), 



(29) 



where P h (k;t) is the linear spectrum at time t, given in terms of the initial spectrum as P h (k;t) = Di(t) 2 P h (k;ti), 
and JC is the function containing the growth functions, which is defined as 



^s(k,r,fi; t) 



1 n 



4 2 i(*) 



r D ch (kr; t)D ch (k y/l + r 2 - 2r^ t) 



1 + r 2 - 2r/i D ch (kr; t)D cb (ky/l + r 2 - 2r>; t) 
fx — r 



r(l + r 2 - 2r M ) D cb (A;r; t)£> cb (A: y/l + r 2 - 2r/x; t) 



(30) 



For the CDM model case (the case without massive neutrinos), i.e. the limit f cb — ► 1, the higher-order growth 
functions become scale-independent as a result of the scale-independence of the linear growth rate. In this case, the 
growth functions are well approximated as 

4ii 2 M 5 



D ch (kr; t)D cb (ky/l + r 2 - 2rfi; t) 7' 



D ch {kr; t)D cb {k^l + r 2 -2rfi; t) 



(31) 



Note that the asymptotic behaviors above are exact only in an Einstein de-Sitter model (r2 m o = 1), but hold an 
excellent approximation for a CDM model at relevant redshifts (5(| |57|. Hence, for the case f cb = 1, Eq. (|2"9"|) recovers 
a well-known expression of the one-loop power spectrum 

p Cdm(M) for the CDM model case (59l. [60j: 

,(22) p(22)„ - k 3 



P 



cb 



prw) 

r-l 



drP^ h (kr;t) 



98(2tt) 2 
d M P c L b (fcVl + ^ 2 ~2r/i; t) 



(3r + 7[i - 10r/z 2 ) 2 



(32) 



.1 {l + r 2 -2r^if 

Thus, for a MDM model, the scale-dependent growth function (Eq. [3D]) has to be solved before obtaining the the 
power spectrum Pj b 2 ' '. Because of this, the exact computation of P^f 1 is computationally expensive. 
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(13) 

Similarly, another one-loop power spectrum, P cb , is formally expressed as 



P 



(13) 



cb, MDM 



(3) 

where the growth function K, s is defined as 



(k;t) 



2k 3 Pt h (k;t) 
(2tt) 2 



dr r 2 P c L b (fcr;i)/Ci 3) (fc,r;t), 



(33) 



K.f(k,r,t) 



D ch (k;t)D ch (kr;t) 2 
2(1 + r 2 ) 



+ 



A ( s 3) (kr, k, kr; t) + ^Af ] (kr, kr, k; t) j 

|sf ) (kr, k, kr; t) + sf ] (kr, kr, k; t) j 
In 



3r 2 

-3 + r 2 (1 - r 2 ) 2 

2 ~ 4r 

-1-r 2 (1-r 2 ) 2 



1 + r 



2r 2 



4r 3 



In 



1-r 
1 + r 



1 -r 



cf\kr,k,kr;t) 



Cf\kr, kr,k;t) 



- 1 [Df (kr, k, kr; t) + Df (kr, kr, k; t) } 



1 + r 2 (1-r 2 ) 2 , 
— In 

2r 2 4r 3 

-l + 3r 2 (1-r 2 ) 2 



1+r 



2r 4 



4r 5 



1-r 
1+1 

In 



Ef\kr,k,kr;t) 



1 - r 



jpf^fcr, fcr,M) 



3r 2 



|pj 3) (fcr, fc, fcr; t) + F 5 (3) (fcr, fcr, fc; t) j 



(34) 



(13) 

Note that P V is evaluated by one-dimensional integration once the linear power spectrum and the growth function 
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are given. For a model without massive neutrinos, the growth functions are approximated as 

1 (A 3 ) r( 3 ) p (3) p(3)\^f 5 1 

D ch (k;t)D ch (kr;tf \ A * ' B& ' C& > D * ^ s J \ 18' 18' 

Therefore, Eq. (|33[) recovers the expression of P^ 13 ) for the CDM model case: 

1 3 poo 

' (13) - ^ 3) ^^2^2W PcLb(M) /o drPkik]t) 

1 + r 



(35) 



P 



cb 



12 



158 + 100r 2 - 42r 4 + — (r 2 - l) d (7r 2 + 2) In 



1-r 



(36) 



Again, an exact computation of P V is computationally expensive. 



Thus evaluations of the one-loop correction spectra, Pgu Z; and P^i! 3 ^, at each k and each time t require high 



dimension integrations, which are numerically time-consuming. Rather, we find that Eqs. (|32[) and (|36[) serve as a 
good approximations to obtain the spectra for a MDM model, if the scale-dependent linear growth rate entering into 
the linear power spectrum P^ h is properly taken into account. We below give the justification. 

First, we study validity of the approximations (|3"Tj) and ([33]) for the nonlinear growth functions in a MDM model. 
Fig. Q] compares the approximations (|3"T]) and (|35[) with the results obtained by numerically solv ing the differential 
equations (|B 1 1) , (|F32|) and (|B4|) that govern time-evolution of the nonlinear growth functions jl03| | . Here we consider 

(2) (3) 

the growth functions ^ (k\,k2) and A s (ki,k2,k^) as representative examples. Also note that we considered /„ = 
0.05, corresponding to the current upper bound, and redshift z = 0, where nonlinear clustering is strongly evolved. 
This figure clearly shows that the fractional errors of the approximations are less than ~ 5% over a wide range of 
wavenumbers we have considered. This agreement implies that the scale-dependence of higher-order growth functions 
are well-captured by the fc-dependence of linear growth rate, D c b(k; t). This level agreement was also found for other 
growth functions, pf } , C| 3) , Df\ Ef \ and P 5 (3) . 

Fig. [5] compares the approximation with the full evaluation of one-loop power spectra, P cb (k) and P cb (k). In 
the left panel, the dashed curves show the result obtained by performing the numerical integrations in Eqs. (|29p 

_____ /2"\ 

and (|33|) where the high-order growth functions arc inserted into the calculations of the Fourier kernels K\ and 
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FIG. 1: The 2nd- and 3rd-order growth functions at redshift z = are plotted as a function of two wavenumbers k\ and &2, 
for a MDM model with j v — 0.05 (m v .to% — 0.6 eV). As representative examples, shown here is the growth functions divided by 
some powers of the linear growth rate: A^P (ki, k2)/[D c t,(ki)D ct>(fe)(5/7)] (left panel), Ag (ki, fe, fci)/[D c b(fci) 2 -D c b(fc2)(5/18)] 
(middle) and A^(ki, k^, k2)/[D c b(ki) 2 D c b(k2)(5/18)] (right), respectively. Note that specific combinations of fci-arguments 
in Af* 1 are chosen because the one- loop power spectrum P c b' ( see ^q. E3) depends on the growth functions of specific 
configurations. The quantities shown become unity for the limit /„ = 0, i.e. a model without massive neutrinos (see Eqs. |3 1 J 
and [3S]), which is shown by the plane in each plot. Therefore the deviations from unity reflect additional scale dependences 
arising from the mode-coupling. It is clear that scale-dependences of the higher-order growth functions are well-captured by 
combinations of the linear growth rate, and the approximations (|3ip and (|35|) hold valid with accuracy better than ~ 5% over 
a range of wavenumbers we have considered. 

K.^, while the solid curves are the results obtained by the approximations (|32[) and ([55)) . Note that we assumed 

f v = 0.05 as in Fig. [T] It is apparent that the absolute values of P^\k) and P^f' ] (k) are slightly overestimated by 
the approximations, because the higher-order growth functions are overestimated by the approximations as implied in 
Fig. [1] The right panel shows the resulting total power spectra of baryon plus CDM perturbations that include up to 
the one-loop corrections: P c b = Pcb + Pcb^ +^ > ib 3 ^ ^ ^ s f° un d that the fractional error of the approximation is smaller 
than ~ 1% on scales up to k — 1/iMpc" 1 , for the case of /„ = 0.05. Thus we can conclude that the approximations to 
compute the one-loop power spectra are sufficiently accurate for our purpose, over ranges of wavenumbers, redshifts 
and neutrino mass scales wc arc interested in. 

Given the results shown in Figs. Q] and we will hereafter employ the approximations (f3"2| and for computing 
the nonlinear power spectrum of total matter at an arbitrary time t. A brief summary is: we first compute the linear 
power spectra of each components at time t, P^ h (k;t), P^(k;i), and P^, v {k\ t), for a desired MDM model, by using 

the publicly available code CAMB [53[. We then compute one- loop power spectra, P^\k;t) and P^f\k;t), using 
Eqs. (|32[) and (|36p . Then, all the spectra are summed up to obtain the nonlinear spectrum of total matter: 

P* h (k; z) = / 2 b [P c L b (fc; z) + F c ( b 22) (k; z) + P^ 3] (k; z)] + 2f ch f„P c L hu (k; z) + f u P^(k; z). (37) 

Before closing this sub-section, we comment on the work of fr3ll |. where a similar method for computing the one- 
loop corrected power spectra for a MDM model was developed based on perturbation theory ignoring the nonlinear 
neutrino perturbations. Although our method is qualitatively equivalent to their method, there are several technical 
differences that may be worth stressing. (6l| first employed the analytical fitting formula for scale-dependent linear 
growth function (also for the transfer function) developed in [54] ]. which is given as a function of neutrino masses 
and cosmological parameters. Then, analytical expressions for the higher-order growth functions and the one-loop 
power spectra were derived. There are several inaccuracies involved in the fitting formula. The formula becomes 
less accurate for small neutrino masses as explicitly pointed out in [62| . More precisely for a case that neutrino(s) 
is massive enough such that the neutrinos become non-relativistic in the radiation dominated era, corresponding to 
m v ^ 0.6 eV (f v ^ 0.05) for a ACDM model, the fitting formula becomes inaccurate because it assumes a continuous 
suppression in the matter growth since neutrinos became the non-relativistic, although the suppression occurs only in 
the matter dominated regime. In addition the fitting formula does not include BAO features in the transfer function. 
Therefore, the use of the fitting formula may under-estimate the ability of future galaxy surveys for constraining 
cosmological parameters, especially dark energy parameters to which the observed scales of BAO peaks are sensitive. 
In contrast, in our method, the nonlinear power spectrum is obtained by inserting the linear power spectrum outputs 
of the numerical Boltzmann solver such as CAMB, which takes into account the scale-dependent growth rate as well 
as BAO features at high precision. Albeit these small differences, [6l| also verified that Eqs. ([3^]) and (f3"o) are good 
approximation. 
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FIG. 2: Left panel: The dashed curves show the one-loop power spectra of CDM plus baryon perturbations, P cb and P cb , 
which are obtained by numerical integrations of Eqs. (|29|l and (|33|l . respectively, while the solid curves show the spectra 
computed using the approximations (|32[l and (|36(l . Note that the y-axis is plotted in the linear scale, and we consider f v = 0.05 
and 2 = 0. For comparison, the rightmost solid curve labelled as "Pcb" shows the linear power spectrum. Right panel: The 
fractional difference of the total matter power spectrum including up to the one-loop corrections is shown in the left panel: 
^ct> L W — Pcb(k) +Pjb 3 ^(^) + -^ib 2 ' W • The approximation is found to be accurate to better than 1% on scales k i$ 1/iMpc -1 . 



III. NONLINEAR POWER SPECTRUM IN A MDM MODEL 



In this section, based on the treatment developed in the previous section, we study effects of finite mass neutrinos 
on the nonlinear power spectrum of total matter. 

In Fig. [3] we show the SPT predictions for nonlinear power spectra divided by the linear power spectra at three 
different redshifts z = 0, 1 and 3, respectively, for a MDM model with f v = 0.01 (m„. to t — 0.12 cV). Nonlinear 
gravitational clustering causes amplitudes of the nonlinear total matter power spectrum to be enhanced, resulting in 
more significant deviations from the linear theory predictions on smaller scales and at lower redshifts. In other words 
our PT model tells the range of wavenumbers and redshifts where the linear theory is valid or equivalently the linear 
theory starts to break down on wavenumbers beyond the applicable range. Our model predictions are also compared 
with the result of an empirical method, which is the halo model approach (hereafter we call 'halofit'). In this model 
the nonlinear power spectrum is obtained by mapping the input linear power spectrum based on the fitting formula 
that is calibrated by numerical simulations for CDM models (63[. Recent studies [l8| employed the halofit method 
to compare the model predictions to the weak lensing measurements for a MDM model, and then derived an upper 
limit on total neutrino mass as TO„ jt ot ~ 0.54 eV (95%C.L.). The halofit power spectra are smaller in amplitudes than 
SPT by up to ~ 10% over a range of scales we consider. Furthermore, SPT more washes out oscillatory BAO features 
than halofit, as pointed out in the previous study J26[, where the SPT results were shown to better reproduce the 
simulations results than the halofit results. 

It should also be noted that SPT eventually ceases to be accurate at smaller scales, and the validity needs to be 
carefully studied by using numerical simulations (e.g., [32| for such a study for a CDM model). Our method may 
be further improved by including the hi gher -order perturbation contributions or using a refined method such as the 
renormalized perturbation theory (e.g. [37j]) or the closure theory method (39| . These are in progress and will be 
presented elsewhere. 

In the left panel of Fig. 0] we compare the two power spectra with and without massive neutrinos, Pf u ^o/ Pf v= o, 
for a fixed O m0 . Note that we show the results for /„ = 0.01 (m^tot = 0.12eV) and /„ = 0.02 (m^tot = 0.24eV) 
to study dependences of the neutrino effect on total neutrino mass, and consider redshift z = 1, the central target 
redshift of WFMOS-like survey. As can be clearly seen, the massive neutrinos imprint characteristic, scale-dependent 
suppression features onto the power spectrum shape. Comparing the linear theory and SPT results manifests that 
the nonlinear power spectrum has increasing suppression on scales k ^ 0.1/iMpc _1 , where the linear theory predicts a 
constant suppression roughly given as P/„^o/-P/„=o ~ — 8/„ [l5j . This enhanced suppression effect can be understood 
as follows. As can be found from Eq. (f32j) and Eq. ([36]). the one- loop power spectra, P^ h and P cb , which give 
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FIG. 3: The nonlinear power spectra for a MDM model with /„ = 0.01. The solid curves show the SPT predictions divided 
by the linear spectra for three redshifts z = 0, 1 and 3, while the dashed curves denote the halofit results. 



nonlinear corrections to the total matter, are roughly proportional to squares of the linear power spectrum, P^ h , and 
therefore the suppression effect on the growth rate is enhanced in the weakly nonlinear regime, compared to the model 
without massive neutrinos. 

The left panel also shows the halofit results. Note that, for this case, the numerator and denominator of P/^o/ Pf v =o 
are both computed by the halofit. Unexpectedly the halofit results fairly well reproduce the suppression features given 
by SPT, although the power spectra themselves show a moderate difference in these two models as implied in Fig. |3] 
The right panel of Fig. 0] shows the results for redshifts z = 0, 1 and 3. The apparent agreement between halofit and 
SPT can be seen only for redshifts z = 1 and 3, and the difference appears clear for for the z = results. Recent 
studies of TV-body simulation in a MDM model also show a similar behaviour of the enhanced neutrino suppression 
[64], [65|. A quantitative comparison among SPT, halofit and TV-body results will be reported elsewhere. 

Is the neutrino effect on total matter power spectrum measurable for a future galaxy survey? To obtain an insight on 
this question, the shaded boxes around the SPT curve with f v = 0.01 display expected 1-er uncertainties in measuring 
band powers of the power spectrum at each wave number bins, assuming survey parameters of WFMOS-like low-z 
survey (see TablcUfor the details). To be more explicit, the fractional errors of measuring the power spectrum, P m (fc), 
averaged over a spherical shell of each radial bin k with bin width Ak are, in an ideal case, given as 



op 



An 2 



V s k 2 Ak 



l + n g P m (fc) 
n g P m (k) 



(38) 



where V s and n g are the comoving survey volume and number density of target galaxies. Note that, for the mea- 
surement errors above, we assumed the Gaussian errors for simplicity, and ignored the non-Gaussian contributions 
(see [66[ for the detailed study). The neutrino suppression appears to be greater than the errors at k ^ 0.06/iMpc -1 . 
Another intriguing consequence of the nonlinear clustering is that the amplified power of P^ L (k) reduces the relative 
importance of the shot noise contamination. Note that in reality matter power spectrum should be replaced with 
galaxy one and we address this issue when forecasting constraints on neutrino masses. Thus, extending available 
range of wavenumber, the constraint on neutrino masses can be improved. 

Finally, it would be worth noting that wiggles in the curves reflect shifts in the BAO peak locations caused 
by the scale-dependent suppression effect of neutrinos. The amount of the modulations, however, is smaller than 
the measurement errors. Hence the uncertainty in neutrino mass is unlikely to largely degrade the power of BAO 
experiments, at least for an expected small f u . 



IV. NONLINEAR GALAXY BIAS 



To model galaxy clustering relevant for actual galaxy surveys we need to further include a galaxy bias effect. 
According to [47} (also see Appendix [Cl for the details), we below describe the modeling of nonlinear galaxy bias in 



12 



1 — Linear ' 

SPT 

Halofit 



E -0.1 



-0.2 




Oh 



£ -0.05 



0.01 



0.1 

k [h Mpc" 1 ] 




o.oi 



o.i 

k [h Mpc" 1 ] 



FIG. 4: Left panel: The fractional difference between mass power spectra with and without massive neutrino contribution. The 
shaded boxes show the expected 1-er errors on the power spectrum measurement for a Stage-Ill type survey of z ~ 1 slice that 
is characterized by the mean number density of galaxies and survey volume, n g = 5 x 10 -4 /i 3 Mpc~ 3 and I4 ur vcy = 1.5 fe~ 3 Gpc 3 
(also see Table U}. The two models of neutrino mass, f v = 0.01 and 0.03 (m v ,tot — 0.12 and 0.36 eV, respectively) are assumed, 
where other cosmological parameters are kept fixed. Right panel: It is shown how the neutrino suppression feature in the power 
spectrum amplitude varies with redshifts, comparing the results for the SPT, linear theory and halofit. 



a MDM model, which is done in a self-consistent manner with the modeling of nonlinear matter clustering presented 
up to the preceding section. 

We assume the local bias: the galaxy distribution at a given spatial position is locally related to the underlying 
matter distribution at the same position, which would be a good approximation at least on large length scales of 
interest. In this modeling the galaxy density fluctuation field, 5 g (x), is given in terms of the matter field S m (x) in a 
Taylor expansion as [67| 

S s (x) = e + Cl S m {x) + ^c 2 8 2 m (x) + ^c 3 S 3 m {x) + ..., (39) 

where c n are the n-th order bias parameters, and e represents the stochasticity of galaxy bias which is a statistical noise 
originating from the fact that the relation between S g and i5 m is not perfectly deterministic. Here the stochasticity is 
assumed to be white noise and be uncorrelated with the density fluctuations, (e^m) = 0. Note that the galaxy bias 
parameters and the stochasticity depend on galaxy type, and vary with time [68l. KM [70j . 

Eq. (|39j) of galaxy bias relation has an analogous form to the perturbative expansion. Therefore we can include 
the bias contribution up to the one-loop corrections in terms of the matter density fluctuations. According to the 
methods developed in [3l|, H3] , the galaxy power spectrum including the one- loop corrections can be computed as 

P g (fc) = b\ [P^L( fc ) + b 2 P h2 . s (k) + b 2 2 P h22 (k)] + N, (40) 

where the functions Ph s ,s and Pb 22 are defined as 

P h2 ,s(k) = 2 J j^I*(q)P&\k-q\)jfHq,k-q), 

P h22 (k) = \ J ^ i*(«)[i*(|fe - g|) - i*(g)]. (41) 

The detailed derivation is shown in Appendix [Cj In Eq. (|40p P^ (fc) is the nonlinear matter power spectrum given 
by Eq. (|37]). and the kernel F s used in Eq. (gTJ) is given by Eq. (fC3"|) . Note that F b 2,5 > and P h22 < at scales of 
interest, 0.01 hMpc^ 1 ^ k 0.2 hMpc^ 1 . Eq. (|4"0|) shows that the nonlinear galaxy power spectrum is modeled by the 
three parameters b\ , b 2 and TV once the matter power spectra are specified for a given cosmological model, where the 
parameters b\ , b 2 and N are redefined from the original parameters in Eq. (|39| and the linear mass power spectrum 
as shown in Appendix [C] (also see [47[ for the detailed derivation) . 

For the limit of very small k, Eq. ([40]) recovers the linear regime result, but with correction term: 

P s (k)^blP^(k)+N. (42) 
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Thus bi acts as an effective linear bias parameter for the pow er sp ectrum and N adds a shot noise contamination arising 
from stochastic bias and nonlinear clustering (also see [Zll, [zl, [zH ) • The terms that depend on 6 2 are proportional 
to the one-loop corrected mass power spectrum give an effect of scale-dependent bias due to the nonlinear clustering. 
These parameters b\ , bi and N change with galaxy type we are working on, so need to be treated as free parameters 
for each galaxy type. In fact, as carefully studied in [3l| . the galaxy power spectrum (|40[) can fairly well reproduce 
the semi-analytic simulation results in the weakly nonlinear regime, if the parameters are properly chosen so as to 
match the simulation results. 

The top panel of Fig. [5] explicitly shows how a scalc-depcndcnt bias in the galaxy power spectrum is modeled by 
62. Note that we consider z = 1 and N = 0, and the spectra plotted are divided by bfP m (k) such that deviation 
from unity represents the effect of scale-dependent bias. Here we consider 62 = ±0.25 and 1.2 as a working example. 
First of all, it is worth noting that the nonlinear galaxy bias of 62 ~ 0.1 — 1 causes a modification in the galaxy 
power spectrum shape at BAO scales, and the effect may need to be taken into account for BAO surveys. For a case 
of 62 < 1, as can be seen from the results of 62 = ±0.25, a positive (negative) 6 2 enhances (suppresses) the power 
spectrum amplitudes increasingly at larger fc, relative to the linear bias case. These features are from the second 
term in the bracket on the r.h.s of Eq. (|40|) because Pb2,s > 0. On the other hand, when 62 > 1, the nonlinear bias 
causes a complicated modification in the spectrum shape, because the third term in Eq. (|40|) becomes dominant over 
the second term at larger k. Note that the third term is always negative, so always suppresses the power spectrum 
amplitudes. 

The results in the top panel imply that, even if a linear bias parameter is well determined, the scale-dependent 
bias may cause a degeneracy with the effect of finite neutrino masses, thereby degrading the ability of future surveys 
for constraining neutrino masses. In particular a negative b 2 causes a suppression in the power spectrum amplitudes, 
similarly to the neutrino effect, so this case may cause a stronger degeneracy. To obtain an insight on this, the middle 
panel of Fig. O studies the neutrino suppression effect on galaxy power spectrum in the presence of nonlinear bias. 
Shown here is the fractional difference of galaxy spectra with and without neutrinos of f y = 0.01, for three cases 
of 62- For comparison, the dashed curves show the results for "matter" power spectra employing linear theory and 
SPT. While the neutrino suppression effect is preserved, the nonlinear bias alters the features in the weakly nonlinear 
regime. This figure shows that a negative (positive) 6 2 weakens (strengthens) the suppression effect. 

The bottom panel shows the dependence on the residual shot noise contamination, given by the term including the 
parameter N in Eq. (|40p . The shot noise term has no wavenumbcr dependence for the power spectrum measurement, 
but the figure implies that the shot noise residual with N = O(10 3 ) (Mpc//i) 3 may significantly alter the power 
spectrum shape over a wide range of wavelengths where the neutrino suppression effect appears. This residual shot 
noise effect arising from nonlinear clustering can be studied by using semi-analytic TV-body simulations where galaxies 
are populated with halos, so may be not so a serious source of systematics in the end (e.g. ff^)- 

Thus nonlinear bias effects cause additional modification on the galaxy power spectrum shape. Therefore uncer- 
tainties in the nonlinear bias parameters need to be properly taken into account in extracting cosmological parameters 
from the measured power spectrum. These will be carefully studied below. 
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FIG. 5: Top panel: The perturbation theory predictions for nonlinear galaxy power spectrum at redshift z = 1, which are 
computed from Eq. (|40[) assuming the three fiducial values of nonlinear bias parameter, 62 = —0.25,0.25 and 1.2, respectively. 
The results are divided by the nonlinear mass power spectrum multiplied by the same linear bias parameter b\ such that the 
deviation from unity represents the nonlinear, scale-dependent bias effect. The positive and negative 62 values, with |&2| < 1, 
cause enhanced and suppressed power spectrum amplitudes on smaller scales compared to the linearly biased power spectrum. 
The model with &2 > 1 causes a complex scale-dependent bias (also see text for the details) . The valid range of linear theory 
and SPT are indicated by the two arrows in the upper horizontal axis (see text for the definition) . Middle panel: The neutrino 
suppression features for the nonlinear galaxy power spectra for different fiducial values of 62- For comparison the two dashed 
curves are the results for mass power spectrum computed from the SPT and linear theory as in Fig. U Bottom, panel: The 
effect of residual shot noise contamination that arises from nonlinear clustering, which is modeled as P g — > P g + N. The three 
solid curves show the results for N = 0, 1000 and 2000, respectively, where bi = 1.51 and hi = 0.25 are kept fixed. 
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V. PARAMETER FORECASTS 

We now estimate the ability of future surveys for constraining neutrino masses when using the SPT model predictions 
to be compared with the measurements. 

A. Fisher Matrix Formalism 

For an actual galaxy redshift survey, there is another nonlinear effect to be taken into account: redshift distortion 
effect due to the peculiar velocities of galaxies. The redshift distortion causes the redshift-space power spectrum to 
be two-dimensional: the galaxy clustering strength is varying as a function of two wavenumbers perpendicular and 
parallel to the line-of-sight direction. The redshift-space power spectrum would be more prominent than the real-space 
one to carry useful cosmological information including dark energy parameters because it contains the geometrical 
distortions in directions both along and perpendicular to the linc-of-sight, the so-called Alcock-Paczynski test [74| 
(also see (7f| [76|). However the distortion effect in the nonlinear regime is not yet fully understood, and a more 
careful study based on high-resolution A-body simulations is needed to develop the accurate modelling (43J. Hence, 
in this paper for simplicity we focus on the one-dimensional real-space power spectrum. This roughly corresponds to 
the monopole power spectrum obtained by averaging the redshift-space power spectrum over the spherical shell of a 
given wavenumber in radius, in combination with the proper weighting as well as with the Finger-of-God compression 
algorithm [48j], as developed in [l3|, [13] ■ Note that, after the spherical shell average, the residual Kaiser's effect of 
redshift distortion [84| behaves like the linear bias parameter. That is, we include only the nonlinear galaxy bias 
effect. 

We can not measure directly the length scale in real space from the observed galaxy distribution; rather we measure 
the angular positions of galaxies on the sky, and the radial position in redshift space. To convert the observed 
position to the real-space position, one needs to assume a reference cosmological model which generally differs from 
the underlying true cosmology. An incorrect mapping causes an apparent distortion in the measured power spectrum, 
known as the geometrical distortion [75l [7r| . Since in this paper we focus on the one-dimensional, real-space power 
spectrum that is given as a function of wavenumber, the wavenumber estimated from the reference cosmology, fc re f, is 
related to the true wavenumber k as 



D v (z) 



ref 



D v (z) 



(43) 



where Dy{z) is the effective distance factor accounting for the spherical shell average in redshift space, and is given 
in terms of the angular diameter distance and the Hubble expansion rate as D v {z) oc \D\(z)j H(z)} 1 ^ (H. The 
quantities with subscript "ref" denote the quantities for the reference cosmology. Further taking into account the 
amplitude shift caused by assuming the reference cosmology, the galaxy power spectrum estimated from a galaxy 
redshift survey, Pj; st , is related to the true spectrum as 



Pf^ref) 



D V (Z) 



P s (k,z). 



(44) 



In order to estimate the accuracies of neutrino mass determination, we adopt the Fisher matrix formalism (e.g. 
see [HI)- The Fisher formalism gives minimal attainable errors on the parameters by means of a set of observables 
considered. However, this method becomes inaccurate in a case that only an upper bound on neutrino masses rather 
than the detection can be obtained for a given survey. In this case we need to take into account the non-Gaussian 
effect of the likelihood, i.e. a sharp cutoff at f u = in parameter space. A more accurate parameter estimation can be 
obtained, e.g. by using a Markov-Chain Monte-Carlo based method [77[. The Fisher matrix formalism is sufficient for 
our purpose, which is to estimate ability of future surveys and to examine the impact of the refined model predictions 
on parameter estimation compared to the linear theory based method. 

The Fisher matrix for the galaxy power spectrum measurement for a given survey is expressed in [78| as 



^galaxy 

a/3 



E 



Vs(Zj) 

4tt 2 



k A dk 



d\nP° st (k;z z ) 91nP K est (fc;zi) 



dpo 



n g (z t )P g cst (fc;Zj) 

n g (z 2 )P„ cst (fc;^) + i 



(45) 



where p a represents a set of free parameters, V s (zi) and fig(zi) are the comoving survey volume and number density 
of galaxies, respectively, at i-th redshift bin defined as [zi — Az/2 : z% + Az/2], and the summation runs over redshift 
slices. Note that P g is given by Eq. (|40|) and the argument k in P^ st is the reference wavenumber k le { in Eq. (I43|) . 
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but wc omitted the subscript for notational simplicity. The partial derivative of the power spectrum with respect 
to parameter p a is computed by infinitesimally varying the parameter p a around the fiducial model assumed, with 
other parameters pp((3 7^ a) being kept to the fiducial values, such that the Fisher matrix estimates the parameter 
accuracies around the fiducial model. 

To compute the Fisher matrix for a given survey we need to specify lower and upper wavenumber bounds in the 
^-integration of Eq. (|45|) . We set k min = 10~ 4 /iMpc -1 , and have checked that choosing the smaller fc m i n little 
changes the results. On the other hand, one should be careful in choosing the maximum wavenumber for each redshift 
slice, fc max (£i), which needs to be chosen from the range of wavenumbers where the model predictions, linear theory 
or perturbation theory, are reliable and accurate. One way to determine k max (zi) is using TV-body simulations in 
comparison with the model predictions. However, high-precision simulations for a MDM model are not yet fully 
explored (see [64|, [65[ for the recent attempts based on the initial pioneer work (zH). Here we simply employ the 
following method for a CDM model in \^ to specify k maK (zi) for each redshift slice: 

P%(q]Zi)dq = C max , (46) 

where is the input linear mass power spectrum at redshift 2^. Since C max is a monotonically increasing function 
with fc max , we will study how a choice of fc max (or equivalently a choice of C max ) affects our results. [321 ] carefully showed 
that, for a CDM model, the SPT results fairly well agree with TV-body simulations up to a maximum wavenumber 
corresponding to C max = 0.18 (0.3) to within up to ~ 1% (3%) accuracy, while the corresponding valid ranges for 
linear theory are given by the smaller values C max = 0.06 (0.13). They also showed that another criterion derived by 
[26l ] seems optimistic, where it was proposed that SPT predictions agree well with simulations up to fc max given by 
A 2 (fc max )= k 3 P(k)/27r\ =ki _ <0.4. 

We also comment that the Gaussian error covariance for galaxy power spectrum is assumed in Eq. (|45p . where the 
power spectra of different wavenumbers are assumed to be independent. Non-linearities of structure formation cause 
correlated errors of different band powers, i.e. non-Gaussian errors for power spectrum measurement, due to the 
nonlinear mode-coupling. The non-Gaussian errors are not negligible even on BAO scales, comparable with neutrino 
free-streaming scales, and are more significant at higher k due to stronger nonlinearities. However at scales of our 
interest the impact of non-Gaussian errors on parameter estimation is expected to be insignificant, so we employ the 
Gaussian error assumption for simplicity (see [66 .. 80] for more detailed studies on the non-Gaussian errors). 

A galaxy survey alone cannot determine all the cosmological parameters simultaneously due to severe parameter 
degeneracies. A useful way to break the parameter degeneracies is combining the galaxy survey constraints with the 
constraints obtained from CMB temperature and polarization anisotropics. In this paper we include information from 
the CMB temperature anisotropy, C; TT , 73- mode polarization, C EE , and their cross correlation, C ; TE , where we use 
the range of multipoles 10 < I < 1500 for C ; TT and C™ and use 2 < I < 1500 for Cf E , respectively. To compute the 
CMB fisher matrix, F™ B , we adopt the noise per pixel and the angular resolution for the Planck experiment that 
were assumed in [82| ■ 

To model the galaxy power spectra and CMB spectra we include all the key parameters that affect the observablcs 
within the CDM and dark energy cosmological framework. Our fiducial model is based on the WMAP 5-year results 
[l|: the density parameters for total matter and baryon are r2 m o(= 0.24), Sl m o/i 2 (= 0.1277), and f2t>o^ 2 (= 0.0223) (note 
that we assume a flat universe); the primordial power spectrum parameters are the spectral tilt, n s (= 1.0), the running 
index, a s (= 0), and the normalization parameter of primordial curvature perturbations, A 2 ^(fco)(= 2.35 x 10~ 9 ) (the 
values in the parentheses denote the fiducial model). We employ the transfer function computed from the CAMB 
code, and note that the primordial spectrum amplitude is normalized at fco = 0.002 Mpc" 1 following the convention 
in [l|. The redshift evolution of dark energy density is given by f2dc(= 1 — f2 m o) and the equation of state parameter 
wq(= —1). When computing the CMB spectra, we further include the Thomson scattering optical depth to the last 
scattering surface, r(= 0.089). For neutrino parameters we assume the standard three neutrino species and vary the 
fiducial value of total neutrino mass, /„. In summary, for a galaxy surveys with N z redshift slices in combination with 
the hypothetical Planck constraints, the model parameters we consider are given as 

Pa = {fl m o,^moh 2 ,^lboh 2 ,w 0l f u ,n s ,a Sl A^ z ,T,bi(z i ),b2(z i ),N(z i )}, (47) 

where z% = z\, 22, %n z - I n total, we include (9 + 37^) free parameters for our Fisher matrix analysis. Note that, for 
the linear theory analysis for the parameter forecasts, we consider (9 + 2N Z ) free parameters (the parameters above 
minus the nonlinear bias parameters 62 The fiducial values of galaxy bias and shot noise parameters change with 
a galaxy survey specification and are described in the next subsection. 

The full Fisher matrix for the joint experiment of galaxy survey and CMB can be obtained simply by adding the 
Fisher matrices: F a p = F^ a * y + F^ B . The unmarginalized error on a given parameter p a is given as u{p a ) = 
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0.2 
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(0.2 < z < 0.4) 






















(0.4 < z < 0.7) 
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0.211 
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0.133 


3.57 




0.65 
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0.263 


0.0194 
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1.61 


1.41 
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0.0177 
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1.31 


Stage-Ill low- z 


1.0 
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4.0 
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0.443 


0.0332 
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1.63 
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0.0124 
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20000 
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1.41 


-0.384 


0.00933 


0.146 


57.2 




1.0 


0.2 


200 


20000 


12.85 


1.51 


-0.345 


0.00594 


0.158 


49.3 


Stage IV 


1.2 


0.2 
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20000 


15.14 


1.63 


-0.299 


0.00383 


0.170 


45.5 


(0.5 < z < 2.1) 


1.4 
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20000 


16.94 


1.74 


-0.242 


0.00217 


0.182 


40.4 




1.6 


0.2 


200 


20000 


18.29 


1.86 


-0.177 


9.96 x 10~ 4 


0.195 


34.8 




1.8 


0.2 


200 


20000 


19.27 


1.99 


-0.096 


2.38 x 10" 4 


0.206 


31.2 




2.0 


0.2 


200 


20000 


19.94 


2.11 


-0.016 


0.06 x 10~ 4 


0.219 


28.4 



TABLE I: Survey parameters that we assume in this paper to make parameter forecasts. The survey parameters are chosen 
such that the surveys fairly well represent the existing survey (SDSS LRG), the near-future planned survey (BOSS), and 
the 5-10 year time-scale future surveys which we call Stage-Ill and -IV surveys, respectively, according to Dark Energy Task 
Force Report [85| . We employ the method described in Appendix [D] (also see text) in order to determine the fiducial values 
of the linear and nonlinear bias parameters bi and &2 for each redshift slice of the respective survey. We also include the 
residual shot noise contamination arising from nonlinear clustering, which is parametrized by N, and we determine the fiducial 
value of each redshift slice according to the method in Appendix [D] The values in the column labelled by k^J^ 3 ^ denote 
the maximum wavenumber up to which the standard perturbation theory (SPT) is expected to be reliable to within a few % 
accuracy compared to JV-body simulation results at each redshift (we determined the fc max values following using Eg. I46p . We 
also show the quantity n g P g (fc max ) at the maximum wavenumber for each redshift slice: if n g P g (fe m ax) > 1, the power spectrum 
measurement is in the sample variance limited regime. 



(Faa) -1 i which corresponds to the accuracy of determining p a when other parameters are perfectly known. On the 
other hand, the marginalize error including uncertainties of other parameters is given as cr(p a ) = [(i r ' _1 ) ceQ ] 1 / 2 , where 
F^ 1 denotes the inverse of the Fisher matrix. The correlation coefficient r between two parameters, p a and pp, is 
also useful to study how the parameters arc degenerate with each other: 

r{p a ,pp)= , {F 1U (48) 



HQ 



If r = +1 (—1), the parameters are totally correlated (anti-correlated), while r = means no correlation between the 
two parameters. 



B. Survey Parameters 



To make meaningful parameter forecasts, we consider survey parameters that fairly well represent future surveys 
being planned or under serious consideration. The hypothetical surveys considered in this paper arc intended to 
resemble BOSS, WFMOS-like survey, and the ideal space-based BAO experiment such as those proposed by JDEM 
and Euclid missions, which are roughly categorized as the Stage-Ill and -IV surveys, respectively, in the DETF report 




The survey parameters are summarized in Table HI The survey area, redshift range and number densities of target 
galaxies were taken from the proposed survey design of each survey. Just briefly, the BOSS-like survey samples 
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Expected marginalized error on total neutrino mass: o-(m„ !t ot) (eV) 







Linear 1% 


Linear 3% 


SPT 1% 


SPT 3% 


+ cr(fi m o) 


= 0.01 


Survey 




C max = 0.06 


Cmax =0.13 


Cmax = 0.18 


Cmax — 0.3 


Cmax — 


0.3 


BOSS 




0.161 


0.111 


0.095 


0.088 


0.082 




Stage-Ill (low- 


z slices alone) 


0.173 


0.123 


0.110 


0.096 


0.082 




Stage-III (low- 


+ high-z) 


0.161 


0.122 


0.107 


0.091 


0.081 




Stage IV 




0.067 


0.059 


0.053 


0.046 


0.046 





TABLE II: Marginalized la error on total neutrino masses, cr(m„ it ot) [eV], expected from each hypothetical survey when 
combined with the Planck and z ~ 0.3 SDSS LRG information. The errors are derived including the galaxy power spectrum 
information over 10 -4 < k < fc max ftMpc" 1 , where fc m ax is determined by Eq. (|46[1 from the input linear mass power spectrum. 
As implied, the linear theory and perturbation theory are expected to be accurate up to the given fc max to within a given 
% accuracy compared to TV-body simulations [32| . For these results we assume m, Ato t = 0.12eV (/„ = 0.01) for the fiducial 
value of neutrino mass, therefore the errors shown roughly correspond to the expected la upper limit on neutrino mass if 
0"(?7V,tot) ^0.12 eV. The last column labelled by "+cr(f2 m o) = 0.01" shows an improvement in the neutrino mass constraint 
for the case C m ax = 0.3 if the prior cr(f2 m o) = 0.01 is added. 

luminous red galaxies (LRGs) over a range of redshifts 0.4 < z < 0.7 extending the SDSS-I and -II surveys. A 
ground-based Stage-Ill survey with optical spectroscopy may be designed to survey galaxies for two different slices: 
one is for galaxies over 0.7 < z < 1.6 with survey area 3200 deg 2 , and the other is for high-redshift Lyman-a emission 
or Lyman break galaxies over 2.5 < z < 3.3. The survey parameters for the Stage- IV type survey are taken from [96[. 
Having multiple rcdshift slices is useful to improve the accuracies of parameter estimation by breaking the parameter 
degeneracies because the sensitivity of each redshift slice to cosmological parameters is slightly different as will be 
shown below (also see [83| for the related discussion). These surveys are complementary to each other in rcdshift 
ranges covered. It is also worth commenting that a high-redshift survey with z > 1 has a potential to explore an early 
dark energy model where dark energy may be more rapidly evolving at higher redshifts than naively expected. 

We further need to specify galaxy bias parameters. However, because we have a limited knowledge on galaxy 
formation, it is difficult to predict galaxy bias parameters with certainty. Here we rather employ a crude method used 
in [23l |75| to estimate the linear galaxy bias parameter b\ for each redshift slice, where 61 is estimated by imposing 
the rms number density fluctuations of galaxies within a sphere of 8/i _1 Mpc radius to be unity: cr 2 R = 1. However the 
LRG bias is relatively well understood based on the existing SDSS sample such as 61 ~ 2.1 in [17| . We assume b\ = 2.1 
for the fiducial value of SDSS LRG bias, from which we compute a correction factor that needs to be multiplied by 
(Tg 8 = 1 to obtain b\ = 2.1 for our fiducial cosmological model. Similarly, for BOSS LRGs, we multiply er 2 8 = 1 
by the same correction factor to estimate the linear bias b\ (see Appendix iDl for more details). The nonlinear bias 
parameter 62 and the residual shot noise parameter TV are more uncertain. We define their fiducial values based on the 
prescription described in Appendix [D] but will study how our results change with different fiducial values of 62 and 
TV. Note that the parameter TV is estimated based on the perturbation theory, but we will employ the same fiducial 
value for the linear theory based forecasts. When TV > l/n g , the residual shot noise contamination is dominant. Our 
survey parameters imply that, for BOSS and Stage-Ill surveys, TV < l/n g . On the other hand, the residual shot noise 
contamination is significant for some redshift slices having higher number densities of galaxies for the Stage-IV survey. 

C. Parameter Forecasts 

1. Summary of constraints on neutrino mass 

Table |TT] summarizes forecasts for the marginalized errors on total neutrino mass for each of hypothetical galaxy 
surveys listed in Table [H combined with the Planck and SDSS LRG information. To derive these errors we determined 
J: max for each redshift slice based on the criteria (|46[) and then included the power spectrum information over 10~ 4 < 
k < k 

max {Zi) /iMpc -1 . We compare the expected constraints obtained when using the linear theory and SPT models, 
over a range of wavenumbers where the respective models seem reliable as indicated from the assumed value of C max 
(the corresponding fc max values when C max = 0.3 for each rcdshift slice arc listed in Table HJ. Note that we assume 
f v = 0.01 (m^.tot — 0.12 eV) for the fiducial value, and the number of free parameters is different in between the 
linear theory and SPT as described around Eq. (|47|) (SPT additionally includes the nonlinear bias parameter 62 for 
each redshift slice). 

It is clear that the use of SPT allows for an improvement in the neutrino mass constraint compared to the linear 
theory results: roughly a factor of 1.3 improvement if SPT can be used up to the maximum wavenumbers where 
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SPT seems reliable to within a few % accuracy corresponding to C max = 0.3. We have checked that the accuracy of 
neutrino mass determination, ^(m^tot), little changes for each survey even if the fiducial value of /„ is varied within 

the current limit f v ^ 0.05. Hence Table [TTI implies the BOSS and Stage-Ill type survey may allow for the accuracy of 
v{ m v,tot) — 0.1 cV, while the Stage-IV survey £r(m„ jta t) ~ 0.05 cV. In particular the expected accuracy for a Stage-IV 
type survey is compatible with the lower limit implied from the normal mass hierarchy. That is, Stage-IV may allow 
for a detection of total neutrino mass at more than l-cr significance; if neutrinos obey the inverted mass hierarchy 
implying the lower limit m v ^ a t ~ 0.1 eV, a 2-er level detection may be achieved. 

Note that the forecasted constraints here are much weaker than those obtained in our previous work [46[. The 
differences are (1) we here consider the one-dimensional power spectrum as the observable rather than the full two- 
dimensional power spectrum in rcdshift space and (2) we include the nonlinear bias parameters. The full analysis 
including the two-dimensional rcdshift power spectrum information will be presented elsewhere (Saito et al. in 
preparation) . 

2. Degeneracy between neutrino mass and other parameters 

To develop a better understanding of the forecasted neutrino mass errors, we study how parameters are degenerate 
with each other and how the degeneracies can be broken when the galaxy power spectrum information ranging from 
the linear to non- linear regime are combined with the CMB information. First, the top panel of Fig. [5] shows the 
unmarginalizcd errors on neutrino mass as a function of the maximum wavenumbcr, for the single z = 1 slice of the 
Stage-Ill low-z survey in Table HI where the linear theory result is compared with the SPT results obtained assuming 
various fiducial values of nonlinear bias parameter 62. For fc max ^ 0.07 ft.Mpc _1 , the neutrino mass constraint does not 
depend on fc max , implying that the constraint is mostly from the CMB information. For the larger fc max the galaxy 
power spectrum is becoming to be more powerful to constrain the neutrino mass due to the increased independent 
Fourier modes. From comparison between the linear theory and SPT results, one can find that the unmarginalizcd 
error on neutrino mass is improved in the weakly nonlinear regime due to the improved signal-to-noise ratio of power 
spectrum measurement, except for the case of 62 < 0. The case of 62 < causes a suppression in the power spectrum 
amplitudes, as implied in the top panel of Fig. [5] As a result the information content of the power spectrum does 
not increase so much in the weakly nonlinear regime compared to the linear theory, although the linear theory breaks 
down in the regime. Thus the neutrino mass constraints are sensitive to galaxy bias parameters or equivalently galaxy 
types. 

The upper-right panel shows the neutrino mass errors marginalized over other parameter uncertainties. Again 
notice that the results are only for one z = 1 slice of the Stage-Ill low-z survey corresponding to the survey volume 
2.1 /i~ 3 Gpc 3 (see below for the full forecast for all the redshift slices combined). Compared to the unmarginalizcd 
errors, the neutrino mass error is significantly degraded due to strong parameter degeneracies. The plot also shows a 
clear plateau feature in the error for /c max ^5 0.1 /iMpc -1 , and then shows a step-like improvement in the error at some 
particular fc max values, which are found to correspond to the BAO peaks. Namely, when the BAO peaks are included 
by increasing fc max , the accuracies of constraining cosmological parameters are dramatically improved by breaking 
the parameter degeneracies via the Alcock-Paczynski test. Comparing the linear theory and SPT results manifests 
that, in contrast to the results for the unmarginalizcd errors, the neutrino mass error does not improve by using 
SPT, due to the significant parameter degeneracies and the addition of nonlinear bias parameter 62- The effect of 62 
can be explicitly studied by the dot-dashed curve, where 62 is kept fixed. Fixing 62 does improve the neutrino mass 
constraints, implying a strong degeneracy between neutrino mass and 62 in the nonlinear power spectrum. However, 
the SPT result with 62 being fixed is still apparently worse than the linear theory extrapolated result in the weakly 
nonlinear regime (although the linear theory breaks down in the regime). This may be understood as follows. As 
discussed, the neutrino mass constraints are sensitive to an inclusion of BAO features which helps break parameter 
degeneracies. However, the nonlinear mode coupling somewhat smooths out BAO features in the weakly nonlinear 
regime, which degrades the constraining power of galaxy surveys in the weakly nonlinear regime. 

The bottom panel of Fig. [6] explicitly studies the Fisher correlation coefficients of neutrino mass with other param- 
eters, r(f v ,p a ), as a function of fc max for the 2 = 1 slice. The neutrino mass appears to be significantly degenerate 
with some parameters such as Q m oh 2 , Wo, fi m o and 62 showing almost perfect degeneracy of \r\ ~ 1. The degeneracies 
show complex behaviors as a function of fc max , where the oscillatory features of r correspond to the BAO features. 

We comment on the parameter N which models the residual shot noise contamination to the power spectrum 
measurement due to nonlinear clustering of galaxies. For the assumed Stage-Ill survey, the residual shot noise 
contamination arising from the nonlinear galaxy clustering is smaller than the standard shot noise l/n s . In addition 
the sample variance gives a dominant contribution to the power spectrum covariance over all the scales we consider, 
k ^ 0.3 ft-Mpc" 1 . Therefore the effect of N is insignificant for the results shown here. However, the genuine effect 
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FIG. 6: Top-left panel: The unmarginalized la error on /„ as a function of the maximum wavenumber fc max , expected from a 
single redshift slice around the centering redshift z c = 1 of Stage-III low-z survey in Table [T] The short-dashed, long-dashed 
and solid curves show how the neutrino mass error changes if the nonlinear bias parameter is changed from the fiducial choice 
&2 = 0.443 to the other choices 62 = —0.443 and 0, where other survey parameters are fixed as shown in Table [I] The neutrino 
mass constraint is sensitive to the underlying &2, i.e. galaxy types. For comparison, the dotted curve shows the linear theory 
result. Note that the two vertical dotted lines show the maximum wavenumbers up to which the perturbation theory is expected 
to be accurate to within the given accuracy for this redshift slice. Top-right panel: The similar plot, but shows the marginalized 
errors on neutrino mass. The step-like features are apparent: the plateau shape is due to strong parameter degeneracies and a 
sudden drop of the error at some particular wavenumbers imply an improvement in the parameter errors because the parameter 
degeneracies can be to some extent broken by including the BAO features with increasing fc max - Compared to the top- left panel, 
the marginalized error is not necessarily more stringent in the weakly nonlinear regime than the linear theory extrapolated error 
due to the stronger parameter degeneracies, although the linear theory ceases to be reliable in the nonlinear regime. Bottom 
panel: The correlation coefficients of neutrino mass with other parameters, r(/„,p Q ), defined by Eq. (|48|l . displaying complex 
degeneracy behaviors as a function of fc max - 



needs to be studied using ./V-body simulations, since this shot noise contamination is not yet fully explored. 

Fig. [7] shows the results combining all the four redshift slices for Stage-III low-z survey (the total survey volume 
V s = 10.24 ft~ 3 Gpc 3 ). The left panel demonstrates the marginalized errors on neutrino mass as in Fig. [71 but as a 
function of C ma x, where fc max for each redshift slice is specified by using Eq. (jifjl) . Compared to Fig. the accuracies 
of neutrino mass determination continue to improve with increasing fc max or adding more galaxy power spectrum 
information; there is no regime dominated by CMB information over the scales we have considered. This is because 
the galaxy power spectra at different redshifts depend on cosmological parameters in different ways, so combining 
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FIG. 7: Left panel: As in the top-right panel of the previous figure, this plot shows the marginalized error on f v obtained 
by combining the four different redshift slices of Stage-Ill low- 2 survey in Table [T] as a function of Cmax, where the maximum 
wavenumber of each redshift slice, fc ma x(2i), is computed using Eq. (|46[1 for an input C max given in the horizontal axis. Note 
that for reference Cmax = 0.1 (1.0) corresponds to fc max = 0.097 (0.266) faMpc -1 at z = 1. The fiducial values of nonlinear bias 
parameters 62 for each redshift slice, which are all positive, are given in Table [I] For comparison, the short- and long-dashed 
curves show the results obtained when the sign of 62 is flipped or assuming 62 = for all the slices, respectively. Right panel: 
The complementarity of different redshift slices is more explicitly studied. The solid curve is same as the solid curve in the left 
panel. The dotted and dashed curves show the results for only one redshift slice with the centering redshifts z c — 1 and 1.2, 
respectively, keeping the survey comoving volume fixed. 

the different redshift information helps break the parameter degeneracies. The solid, short- and long-dashed curves 
compare the results for different fiducial values of 62: In the first case we adopt the fiducial values of 62 given in 
Table [l] while in the second case we multiply the minus sign in the fiducial values of 62 in Table U (all the 62 values 
are negative for this case). In the third one, we set 62 = for all redshift slices. As also discussed in Fig. the 
neutrino mass determination accuracies are found to be sensitive to the fiducial values of 62 , or equivalently to galaxy 
types targeted for future surveys. However the differences due to different values of 62 become milder by combining 
different slices. 

The SPT results can be compared with the linear theory result (the dotted curve). As in Fig. the nonlinear 
regime suffers from severe parameter degeneracies, yielding less stringent parameter constraints than naively expected 
by linear theory. For this reason the parameter forecasts in the previous studies may be somewhat too optimistic, if 
the forecasts are derived based on the linear theory and the linear bias parameter (e.g., [H, HH EH )■ However, we 
again note that the full information on galaxy clustering is inherent in the two-dimensional redshift space, while we 
consider the one-dimensional power spectrum in this paper. 

The usefulness of combining different redshift slices is explicitly shown in the right panel of Fig.[7J The plot compares 
between the results of different redshift slicing, where the survey volume is kept fixed to V s ~ 10.24 h~ 3 Gpc 3 . However, 
note that the maximum wavenumber fc max is different for different redshift slices, therefore the effective survey volume 
is different. The solid curve is the result of our fiducial Stage-Ill low-z survey, while the dotted and dashed curves 
are the results assuming single redshift slice which have different centering redshifts z c — 1.0 and 1.2 with width 
Az = 0.2, respectively. The single redshift slice cases correspond to survey areas f2 s = 15900 and 13500 deg 2 , 
respectively, compared to the fiducial area f2 s = 3200 deg 2 over redshift range 0.7 < z < 1.6. It is clear that the 
neutrino mass constraint is improved by combining the different redshift slices. Also, comparing the dotted and 
dashed curves clarifies that a choice of redshift slices affects the constraining power in the weakly nonlinear regime. 

3. Impact of massive neutrinos on dark energy constraints 

The primary science goal of future surveys is constraining the nature of dark energy via the BAO experiment. 
However, the dark energy constraints may be biased if the model fitting ignores neutrino mass contribution. Fig. [H] 
presents the marginalized error ellipses in a sub-space of the neutrino mass /„ and the dark energy equation state 
parameter wq for the Stage-Ill and -IV surveys, respectively. Note that the dark energy constraints shown here are 
from both the BAO peak locations and the power spectrum amplitude information. There appears to be a significant 
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FIG. 8: Forecasted lo error ellipses in (wo, /„)-subspace for Stage-Ill (top two panels) and Stage IV (bottom), respectively. 
The outermost, intermediate and innermost contours show the results assuming C m ax = 0.13, 0.3 and 0.7, respectively, that 
correspond to the higher fc max in each redshift slice. Again note that C ma x = 0.13 and 0.3 roughly correspond to the maximum 
wavenumbers that the linear theory and SPT are reliable with a few % accuracy. The case of C max = 0.7 may be feasible if the 
refined model of nonlinear power spectrum can be used (see text for the details). The dashed contours show the errors when 
the prior of o(b2) = 0.01 (left panels) or cr(fi m o) = 0.01 (right panels) is added. 



correlation between wq and neutrino mass as expected. For example, a model with wq > — 1 or greater fideo yields 
smaller amplitudes in the galaxy power spectrum, because such a model causes dark energy to be more significant 
from earlier epochs and therefore the greater cosmic acceleration suppresses the clustering growth rate for the CMB 
normalization of linear power spectrum amplitude. This dark energy effect can be compensated by lowering the 
neutrino mass (i.e. the smaller /„) that leads to a less suppression in the power spectrum amplitudes at the larger 
k. One can also find that having larger fc max (equivalently larger C max ) yields more stringent constraints on these 
parameters. In particular, it should be noted that a Stage-IV type survey may allow for a stringent test of neutrino 
mass, even from the ID power spectrum information over a range of wavenumbers where SPT seems reliable. 

In Fig. [8] we also study how the parameter constraints are improved by adding an external prior of 62 or O m o. 
These priors may be delivered from the galaxy bispectrum analysis [8lj . the SN survey and/or weak lensing surveys 
[80l ]. Adding the priors shrinks areas of the error ellipses, because f2 m o and 62 are degenerate with neutrino mass and 
dark energy parameters in the galaxy power spectrum as implied in Fig. [5] In particular, for a Stage-Ill type survey, 
the prior of precision a(fl m o) ~ 0.01 can efficiently break the TO„,tot-wo degeneracy, thereby yielding the accuracies 
of <j(m Vyt ot) — 0.1 eV and a(wo) ~ 0.05, respectively. For a Stage-IV type survey, the constraining power is already 
sufficient, so such a prior does not much help improve the parameter constraints. 

A more important question is how the uncertainty of neutrino mass affects dark energy constraints from future 
galaxy surveys. Table Hill addresses this issue. First, comparing between the third and fourth columns clarifies that the 
accuracy of wq determination is affected by including neutrino mass parameter in the model fitting. If neutrino mass 
is ignored, the error of wq is apparently tightened by a factor of 1.2-1 A for the galaxy surveys we consider here. It 
should be noted that the tighter constraints correspond to a case that the neutrino mass is sufficiently well determined 
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Survey 


range of k 

(Cm ax) 


a (wo) 


a(w ) 
ignoring /„ 


5 u>o 
ignoring /„ 


a(w ) 
+a(Q m0 ) = 0.01 


5 wo 

+a(Q m0 ) = 0.01 


BOSS 


Linear 3% (0.13) 


0.1522 


0.0978 


0.0090 


0.0507 


0.0262 




SPT 3% (0.30) 


0.0768 


0.0603 


0.0141 


0.0435 


0.0243 


Stage-Ill 


Linear 3% (0.13) 


0.1935 


0.1067 


0.0060 


0.0503 


0.0255 




SPT 3% (0.30) 


0.1103 


0.0801 


0.0125 


0.0476 


0.0254 


Stage-IV 


Linear 3% (0.13) 


0.0398 


0.0375 


0.0113 


0.0311 


0.0176 




SPT 3% (0.30) 


0.0245 


0.0223 


0.0206 


0.0226 


0.0223 



TABLE III: The impact of massive neutrinos on determination of dark energy equation of state parameter, wo- Note that 
f u = 0.01 is assumed for the fiducial model. The fourth and fifth columns, labelled as "ignoring /„", show apparently tighter 
constraints and biased best-fit values of wo caused when galaxy power spectrum models without neutrino mass parameter is 
fitted to the true spectrum with /„ = 0.01, respectively, for each galaxy surveys. The sixth and seventh columns, labeled as 
"+cr(f2 m o) = 0.01", show the similar results when adding the prior <x(O m o) = 0.01. 

by a laboratory experiment (i.e. in this case f v is no longer a free parameter in the galaxy power spectrum). 

More importantly, ignoring neutrino mass in the model galaxy power spectra likely results in a biased best-fit value 
of Wq. According to the method described in Appendix [El the column labelled as "Swq ignoring /„" estimates an 
amount of the possible bias, that is, the difference between the input wq and the best-fit value obtained from the model 
fitting without neutrino mass parameter: WQ est ~ Rt = —1 + Swq. To be more explicit, we here estimate the bias caused 
when the template of galaxy power spectrum assuming /„ = is fitted to the observed spectrum having the true 
neutrino contribution of f„ = 0.01 (m u ,tat = 0.12 eV). The table shows a positive bias 5wq: iz;Q est- * > WQ nput = — 1, 
because a model with wq > — 1 predicts galaxy spectra with smaller amplitudes due to the suppressed growth rate, 
which mimics the neutrino suppression effect inherent in the (presumably here) measured spectrum. For BOSS and 
Stage-Ill type surveys, the bias is not significant because \5wq/<j(wq)\ < 1, while a Stage-IV type survey may suffer 
from a significant bias as \Swq\ ~ <r(wo). Table Hill also shows that a 1%-level prior of f2 m o helps reduce the statistical 
error a(wo), but also make the systematic bias more significant at the same time. 

Fig. [9] more nicely illustrates the impact of neutrino mass uncertainty on dark energy parameter estimation from 
future galaxy surveys, showing the projected error ellipses in (w , ^doo) _ pl ane - It is clear that ignoring /„ leads to 
model fitting apparently with smaller error ellipses and biased best-fit values for these parameters. In particular, 
for a Stage-IV type survey, the biased best-fit dark energy model confined by the 1-er statistical error bounds may 
happen to be outside from the underlying true model (the input value wq = — 1 in our case). The amount of bias 
would become greater for the greater values of true neutrino mass. Thus Table IIIII and Fig. [9] imply that neutrino 
mass contribution is not negligible and needs to be included in the model interpretation for future galaxy surveys in 
order not to have too optimistic and biased dark energy constraints. 

Note that the parameter biases studied here are mainly from the power spectrum amplitude information. If the 
dark energy parameters are estimated from BAO peak locations being marginalized over a sufficient number of 
nuisance parameters that include power spectrum amplitude parameters (3fjj . the dark energy parameter biases can 
be minimized, although the constraining power is significantly weakened. This is beyond the scope of this paper, but 
would be worth carefully studying. 

VI. SUMMARY AND DISCUSSION 

Following our earlier work (4o| . in this paper we have developed a method for computing nonlinear power spectra 
of total matter and galaxies in a mixed dark matter (MDM) model (a model with CDM plus finite-mass neutrinos) 
based on standard perturbation theory (SPT) approach. In particular we have carefully examined the validity of 
approximations employed in our approach. 

For our fiducial approach, we include only the linear-order neutrino perturbations to compute the nonlinear power 
spectrum, where nonlinear clustering is driven by the nonlinear growth of CDM plus baryon perturbations. Our 
approach is motivated by the fact that the neutrino free-streaming scale is sufficiently large for small neutrino mass 
scales consistent with the current limit (m„ ito t ~ 0.6 eV) and the neutrinos are expected to more stay in the linear 
regime than CDM plus baryon. We carefully studied the validity of this assumption as briefly summarized in the 
following (see Appendix . By solving the hierarchical Boltzmann equations of neutrino perturbations including 
the nonlinear gravitational potential contribution (due to the nonlinear CDM and baryon density perturbations), we 
indeed found that the amplitudes of neutrino density perturbations are enhanced by up to 10% at the weakly nonlinear 
scales up to k ~ 0.5Mpc _1 . Thus, although the neutrino perturbation is indeed affected by nonlinear clustering, the 
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FIG. 9: The projected la error ellipses in (wq, f2DE)-plane for Stage-Ill (left panel) and Stage-IV surveys (right panel), 
respectively. Note that C ma x = 0.3 is assumed. The solid contours in each panel shows the result for our fiducial method 
where neutrino mass contribution to the galaxy power spectrum is properly taken into account and the errors on dark energy 
parameters are derived by marginalizing over other parameter uncertainties. The dashed contours and the triangle or square 
symbols show the worst-case results: apparently tighter constraints (smaller error ellipses) and biased best-fit values may be 
caused if neutrino mass contribution is ignored in the model galaxy power spectra. The triangle and square symbols show the 
biased values when the underlying true cosmology has f v — 0.01 (m„ jt ot = 0.12 eV) and 0.02 (0.24 eV), respectively. For a 
Stage-IV type survey, ignoring neutrino mass may cause a false best-fit model that is away from the true model wo — — 1 by 
more than the 1-a statistical errors. 

contribution to total matter power spectrum can be safely ignored to less than a sub-percent level, for neutrino mass 
scales of interest, because the neutrino perturbation contribution to total matter clustering is suppressed by another 
small factor /„ = f2„o/f2 m o whose current limit /„ ^ 0.05. 

Further we carefully studied the higher-order growth functions of CDM plus baryon perturbations. Compared to the 
CDM case, the finite-mass neutrinos cause a scale-dependent suppression in the clustering growth rate, and therefore 
the higher-order growth rates generally have complicated scale dependence. That is, the additional nonlinear mode- 
coupling between perturbations of different wavenumbers arise via the growth rates in a MDM model. We numerically 
solved the differential equations of the higher-order growth rates, and found that the higher-order growth rates are 
well approximated by the power of the linear growth rate (see Fig. [TJ). 

As a result the nonlinear power spectra can be approximately given by rather simple forms (see Eqs.|32j. |36| 
and [37]) similarly as in the SPT approach for a CDM model. The equation (|3~T|) is very useful in a sense that the 
nonlinear power spectrum at a given redshift z can be computed from the linear transfer functions of CDM, baryon 
and neutrino perturbations at the redshift z, which are standard outputs of the publicly available codes, CMBFAST 
or CAMB. As in (46[, we found that the neutrino suppression effect on the total matter power spectrum amplitude 
is enhanced in the weakly nonlinear regime than in the linear regime (see Fig. @|. Note that the empirical halofit 
approach shows 10%-level deviations from the SPT results in the weakly nonlinear regime, although it qualitatively 
captures the neutrino effect in the nonlinear regime (see FigUJ. 

Thus we believe that our approach gives reliable, accurate model predictions for the nonlinear matter power spec- 
trum in a MDM model over a wider range of scales, where the perturbation theory is valid, than the linear theory. 
Also important is the SPT approach can explicitly tell the scales and rcdshifts where the linear theory ceases to be 
accurate or breaks down. However, simulation based studies are definitely needed to test and/or calibrate the SPT 
predictions. An iV-body simulations for a MDM model is still challenging, but encouragingly the initial attempts are 
being explored in [64|, ■ Alternative approach to refine the analytical modeling is to extend the SPT approach by 
including higher-order loop corrections. Recently there are several efforts made in this direction for a CDM model: 
the Time-RG formalism [4(J , the renormalized perturbation theory (RPT) (3?j] and the closure theory approach (39j 
some of which show remarkable agreement with TV-body simulations over a wide range of B AO scales [12, EH ■ For 
example, an attempt to extend the PT approach for a MDM model has recently been made in the framework of the 
Time-RG formalism [§(J. Given the similarity between our approach and the PT of CDM model as described above, 
we hope that our method can be straightforwardly extended to the improved nonlinear modeling. This is our future 
project, and will be presented elsewhere. 
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Another interesting result of this paper is we developed a method to compute the nonlinear galaxy power spectrum 
in a MDM model by taking into account the nonlinear biasing effect in a self-consistent manner within the SPT 
framework. As given by Eq. (|40|) . the nonlinear galaxy power spectrum is modeled, in addition to cosmological 
parameters, by introducing the linear and nonlinear bias parameters, b\ and 62, and one additional parameter to 
model the residual shot noise contamination N . Once again, although the validity of SPT approach needs to be 
tested by simulations, our SPT approach is built on the physical foundation of large-scale structure formation and 
therefore expected to be reasonably accurate in the weakly nonlinear regime where the SPT approximately works out. 

After formulating the nonlinear galaxy power spectrum, we then estimated the ability of future galaxy surveys for 
constraining neutrino masses from the power spectrum information over scales ranging from the linear regime to the 
weakly nonlinear regime. In this paper for simplicity we focused on the real-space power spectrum, i.e. ignored the 
redshift distortion effect, because the nonlinear distortion effect, the Finger-of-God effect, is not yet fully understood 
even in the weakly nonlinear regime. We found that the accuracy of neutrino mass constraint is indeed improved by 
including the power spectrum information up to the weakly nonlinear regime compared to the linear regime, by a 
factor 1.3, for all the planned BAO surveys (see Table Hl|) . However, the improvement is not so significant because 
of severe parameter degeneracies in the nonlinear regime (see Figs. |6]and[7|)- Thus the neutrino mass forecasts in 
the previous studies may be too optimistic if the forecasts are derived assuming the linear bias and the linear theory 
modeling. Nevertheless it should be noticed that Stage-Ill and -IV type surveys may allow for the neutrino mass 
constraints to accuracies of ~ 0.1 and 0.05 eV, respectively, even from the ID power spectrum information. 

We also studied how the finite-mass neutrinos affect the ability of future surveys for constraining dark energy 
parameters. A change of dark energy parameters such as wq > — 1 from cosmological constant model also causes a 
suppression in the galaxy power spectrum amplitudes, because the growth rate of mass clustering slows down due to 
the greater cosmic accelerating expansion. Thus the dark energy constraints are likely correlated with neutrino mass 
in the galaxy power spectrum (see Fig. |5J), if the power spectrum amplitude information is included in parameter 
estimation. In particular we pointed out that, if neutrino mass parameter is ignored in the model fitting, the best-fit 
dark energy parameters can be biased. For a Stage-IV type survey, the bias may be greater than the statistical 
uncertainty: a false evidence of wq 7^ —1 may be implied by the neutrino mass uncertainty, even if the true model 
has wq = — 1. Thus our results suggest that the neutrino mass contribution needs to be taken into account for future 
BAO surveys and to be marginalized over in order to obtain an unbiased constraint on dark energy parameters. 

We believe that the SPT modeling of galaxy power spectrum can be a more physically motivated model than other 
empirical approaches such as the halo model approach or the method where nuisance parameters such as Qnl in [13 
were empirically introduced to model the nonlinear effects including the nonlinear bias effect. The method developed 
in this paper allows us to model the nonlinear galaxy power spectrum self-consistently within SPT formulation without 
introducing empirical nuisance parameters. Hence we hope that the use of SPT model allows an unbiased extraction 
of cosmological parameters from the measured galaxy power spectrum by marginalizing over the bias parameters, as 
long as the analysis is restricted to scales where SPT is valid. We are planning to apply our method to the SDSS LRG 
power spectrum. For the SDSS power spectrum measurement done in [ijj , the redshift distortion effect is supposed to 
be removed by using the Finger-of-Got compression algorithm (4^|. Note that the residual Kaiser's effect of redshift 
distortion is absorbed in the linear bias parameter after the spherical shell average of galaxy power spectrum in 
redshift space. Therefore the LRG power spectra are appropriate to compare with the SPT model predictions studied 
in this paper. We will address how the use of our SPT model changes the neutrino mass constraints as well as other 
cosmological parameter estimation as a function of the maximum number fc max , compared to the previous results. 
This is now in progress and will be presented elsewhere. 

There are several other applications of our method. First is gravitational lensing effects on CMB or distant galaxy 
images, which arc sensitive to total matter distribution and therefore known as a powerful probe of neutrino mass 
being free of galaxy bias uncertainty (e.g. [13, HH). These lensing signals are affected by nonlinear clustering, but 
the effect for a MDM model has not been fully explored. Secondly, the formulation developed in this paper tan 
be straightforwardly extended to studying the higher-order correlations of total matter and/or galaxy distribution, 
based on the SPT approach. The higher-order correlations are expected to be very powerful to improve cosmological 
constraints when combined with power spectrum information, and especially to break degeneracies with galaxy bias 
parameters for a galaxy clustering case. 
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Throughout the paper, we assumed that neutrino perturbations stay at linear level, and contribute to the higher- 
order CDM plus baryon perturbations only via the effect on the growth rate. We then simply used the result of 
neutrino perturbation in linear theory, 8^. This assumption is essential for our formalism. In this Appendix, we 
discuss the validity of this assumption in some details. 

Rigorously speaking, the higher-order Boltzmann equations for massive neutrinos must be solved for a quantitative 
estimate of the nonlinear effect on neutrino perturbations. However, there are at least two important facts that 
simplify the analysis. One is that the nonlinear gravitational instability is mainly driven by the CDM plus baryon 
perturbations, which have a dominant contribution to the total matter density, / c b ^ 0.95. Another important fact is 
the presence of the neutrinos' large free-streaming, which prevents the neutrinos from clustering together with CDM 
plus baryon on scales smaller than the neutrino free-streaming scale. Hence, to a good approximation, the impact 
of nonlinear clustering on the neutrino perturbations can be estimated from the nonlinear gravitational potential <fi 
driven by the nonlinear CDM plus baryon perturbations, just ignoring the higher-order neutrino perturbations. Note 



where P\ and P b describe the nonlinear CDM plus baryon density perturbations and are calculated from Eqs. (j32|) 
and ()36|) . Provided P c b for a given cosmological model, we numerically solve the linearized Boltzmann hierarchies, 
Eqs. (f9jl- (fTTj) coupled with Eq. (|Alj) . and obtain the solutions for at a given redshift. We have used the CAMB 
code to implement this approach, modifying the corresponding parts in the code. Note that in the Poisson equation 
given above, the nonlinear corrections to the power spectrum P c b are calculated assuming the linearity of neutrino 
perturbations. In this respect, our approach is not self-consistent, but is sufficient for our purpose to estimate the 
impact of the nonlinear clustering. In fact, the effect is found to be sufficiently small for scales of interest as shown 
below. Furthermore, if necessary, the correction to the CDM plus baryon perturbations due to the nonlinear neutrino 
perturbations can be computed iteratively in a perturbative manner. 

Fig. [10] shows the fractional difference between the linear-order neutrino density perturbation 6^ and the nonlinear 
perturbation obtained from the treatment mentioned above, where (5^ L is calculated by inserting the solution 
for hito Eq. (fT4"|) . We here chose a rather large neutrino mass, /„ = 0.05, close to the current upper bound. The 
plot clearly shows that nonlinear gravitational potential indeed enhances the neutrino perturbation by up to ~ 10% 
on scales where the PT is presumed to be valid. Since the contribution of neutrino perturbation to the total power 
spectrum always involves small additional factor /„ (see Eg. I17[) , the result implies that influence of non-linearity on 
the total matter power spectrum is much more reduced. As a result, we found that the amplitude of P m (/c) increases 
only by 0.01% compared to that obtained using the method of our paper. This effect gets even smaller as < 0.01% 
when /„ < 0.05. Hence, the error caused by the assumption that the neutrino perturbations stay at linear level is 
safely negligible, compared to the measurement errors at a percent level for a future survey. 

Finally, we briefly comment on the recent work in [92| . They discuss the effect of higher-order neutrino perturbations 
just treating the neutrinos as fluids with pressure. Strictly speaking, neutrinos cannot be treated as fluid, and the 
higher-order effect of moment hierarchy should be taken into account in a self-consistent way. Moreover, their 
formulation heavily relies on the assumption that neutrino perturbations stay at the same order as in the case of 
CDM plus baryon fluctuations, which is manifestly violated in the presence of the neutrino free-streaming. Even at 
the linear-order level, 3> 5v at the scales smaller than the neutrino free-streaming. Nevertheless, their results 
are qualitatively similar, and agree well with those examined here. 



APPENDIX A: NONLINEAR EFFECT ON NEUTRINO PERTURBATIONS 
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APPENDIX B: HIGHER-ORDER GROWTH FUNCTIONS IN A MDM MODEL 



In this appendix, we summarize the basic equations for higher-order growth functions in a MDM model defined in 
Sec III C[ which were used for the analysis presented in Sec. Ill Dl 
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FIG. 10: The fractional difference between the linear-order neutrino density perturbations 5„ L computed using the approach 
described in text. We assume f v = 0.05 that corresponds roughly to the upper bound of current observational upper-bound. . 



Let us consider the second-order growth functions, and B\ Z) , defined in Eq. (j26|) . From the perturbation 
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equation for second-order quantity S ch (see Eq. (|25|) ). the governing equations for A s and B s are obtained, and we 



have 
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Note that in numerically solving the above equations, we retrieve only the inhomogeneous part of solutions so that 
the solution consistently approaches zero when going back to an initial time t — > 0. This treatment just corresponds 
to picking up the growing-mode solution consistently, since the source terms of the evolution equations involve the 
growing- mode solution of linear perturbations. 

Next write down the governing equations for third-order growth functions defined in Eq. (f2"T|) . I s (fci, k 2 , fc 3 ) (X = 

A, B,C, D, E, F), shortly abbreviated as Ij . To do this, we first derive the perturbation equation for third-order 

quantity S ch . From Eq.©, substitution of the linear and second-order solutions <5^ b ' and ci J leads to 
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where the quantities, ai j23 , a 23 i and /?i,23 respectively indicate fc 2 +fc3), a(k 2 +k 3 , k±) and /3(fei, fc 2 +fc 3 ). Then, 
comparing the formal solution (|27[) with the above equation, we obtain the evolution equations for the third-order 

growth functions ij 3 ^ : 
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Again, the above equations must be solved just retrieving the inhomogeneous part of the solution. Here, the source 
functions, Sf 2 3 , {I = A, B,C, D, E, F), are the scale- and time-dependent functions consisting of the linear and 
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second-order growth functions. They are given by 
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In the above, the functions and i?g are the second-order growth functions which appear in the solution of 

(2) (2) (2) 

second-order velocity divergence 9^. h ' . These functions are related to the functions A g ' and B y g through 



Af\kx,k 2 ;t) = u±y f^ L > Dch (k 2 ; t) - H^Af (k u k 2 ;t), B { f\k u k 2 ;t) = -H~ l Bf k 2 ; t). (B6) 
d In a 

Finally, we note that in the limit of / c b — > 1 (i.e., case of massless neutrinos), there exist no free-streaming scales, 
and the linear growth function D c b becomes independent of scales. From Eqs. ()Bll) . (|B2I) and (|B4|) . this readily implies 
that all the second- and third-order growth functions become scale-independent. Then, employing the Einstcin-dc 
Sitter approximation, the analytical expressions for higher-order growth functions can be systematically obtained. In 
the Einstein-de Sitter approximation, all the calculations done in the Einstein-de Sitter universe are extended to apply 
to the other cosmological model by simply replacing the linear growth function in the Einstein-de Sitter universe with 
the one in the underlying cosmology. The detailed discussion on the validity of the Einstein-de Sitter approximation 
is given in (56|, H3, \EM ■ 

As a result, higher-order growth functions in the / c b — > 1 limit arc analytically expressed as 
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for the second-order growth functions, and 

\ A f\Bf\cf\Df\Ef\F^\^{^,l-- 1 -- l --l-,-^-\ D x {tf. (B8) 
\ 8 » s ' s > 5 ' s > s / \l8'l8 6 9 21 63 J w y J 

for the third-order growth functions. For reference, in Fig. [TT]we show nonlinear growth functions as in Fig.[TJ 



APPENDIX C: REPARAMETRIZATION OF BIASING PARAMETERS 



In this appendix, we review the reparamctrized biasing parameters proposed by Ref. [47| . In this treatment, the 
galaxy power spectrum can be consistently related to the matter power spectrum calculated from SPT. 

The starting point is that the fluctuation of galaxies is expanded in Taylor series assuming the local biasing 
prescription. In the local biasing scheme, the galaxy density field at a given position is described as the local function 
of matter fluctuation at the same position. In Fourier space, the galaxy density field is described as 



5 g (k) 



ci5 m (k) + y J 5 m (q)5 m (k - q) 

c 3 f d 3 q 1 d 3 q 2 



6 



(2tt) 6 



M<h )M9 2 )Mfc - Qi - Q2) + e(fc) + O(Sif ), 



(CI) 



where the Cj's are the biasing parameters. The quantity e represents the residual random field which cannot be 
represented by the matter fluctuations. We assume that randomness of e is described by a white noise, and is 
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uncorrelated with <5 m , i.e., (e 2 ) = Nq and (e5 m ) = 0. Then, the galaxy density power spectrum P g up to the one-loop 
level is calculated as 

P g = c 2 P^(k)+Lc 3 a 2 + ^c lC2 a 2 )p^k) 



+2c lC2 J P%(q)I*(\k-q\)jf\q,k-q) 



(27T) 



(2tt) s 

P^(q)P^(\k-q\)+N , (C2) 



d 3 q 



where the constant parameter a 2 is defined as a 2 = J d 3 q P^(q) / (2t:) 3 , and the function Tg 2 \k,k') is the Fourier 
kernel of the second-order density perturbation: 



5 1 k ' k f k k \ 2 ( k • k 



While the calculation in the above is exact up to the fourth-order in density, due to the truncation at finite order, 
the expression (|CJ2|) suffers from several unphysical behaviors such as an apparent divergence and anomalous low-fc 
power [71j . To remedy this, Ref. [47j proposed a way to regularize the expression (|C2[) by reorganizing several terms 
and reparametrizing the biasing parameters. In this treatment, the first line of Eq. (|C2|) is rewritten as 



Cl c 3 a 2 + f iCl c 2 a 2 ) Pi 



2 2 ^8 r 



P^(k). (C4) 



Further, the apparent divergence arising from the third line is absorbed by redefining the parameter N as 



Then, the galaxy power spectrum is re-expressed as follows: 

P g (fc) = b\ [P^ L (k) + b 2 P h2 . s (k) + b 2 P h22 (k)] + N, (C6) 

iW(*0 = 2 / ^(9)^(|fe-g|)^i 2) (g,fe-9), (C7) 

P h22 (k) = \ J i*(«)[i*(|fc - 9 |) - i*(g)]. (C8) 

As a result, the galaxy power spectrum in the weakly nonlinear regime can be described with the only three parameters, 
bi,b 2 and N. Recently, the validity of the expression (|C8[) has been examined in some details in Ref. [3l|. They 
reported that this reparametrization scheme can fit well to the power spectrum of halos and galaxies in millennium 
simulations, and the cosmological parameters can be correctly estimated using Eq. (|C8j) as a template in an unbiased 
fashion. 



APPENDIX D: NONLINEAR BIAS PARAMETERS BASED ON HALO-MODEL 

We here summarize how to determine the fiducial values of the biasing parameters b\ , b 2 and N listed in Table [H 
which are used in the Fisher matrix analysis in Sec. IV CI 

Following the treatment in Refs. [HI, [75[, we determine the linear biasing parameter bi at a given rcdshift so that 
the condition ag tg (z) — 1 is satisfied, where we define 

a 8 , g (z) = 6i(^( ? )^7^ + ^f-, (Dl) 

with F m (z) = —din Di(z) / d\n(l + z). The function Di(z) is the linear growth rate for ACDM model, which we 
compute from the fiducial cosmological parameters just setting f v = 0. For SDSS LRG and BOSS surveys, their 
target samples are LRGs whose clustering properties arc relatively known from the observations, and the linear 
biasing parameter is measured as b± ~ 2.10 at z = 0.3. Hence, when considering these surveys, we simply adopt this 
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value, and the linear biasing parameters at different rcdshifts arc determined from (|D1|) just rescaling the condition 
°8,g( z ) = 1 to a s , g {z) = er 8ig (0.3) with &!(0.3) = 2.10. 

The non-linear biasing parameter b 2 in the expression (|C8|) is related to the original parameters Ci in Eq. (|Cl| as 



b 2 = (D2) 

ci 

The biasing parameters C\ and c 2 can be estimated from the halo-model approach (e.g., [93|). According to this 
prescription, we obtain 

1 f°° 

Ci = — dM n b (M,z)b}(M,z)(N) M , (D3) 

n S JM min 

where the function rih(M, z) is the halo mass function for the given mass M and redshift z, and the quantity b^(M, z) is 
the halo biasing parameter. The expectation value (N)m is the so-called halo-occupation distribution, which describes 
the mean number of galaxies per halo with mass M. We here set (N)m = 1 for simplicity. We adopt the Sheth and 
Tormen formula for mass function n^{M, z) (94j : 

n h (M, z) = -E^^fiu) ; /(„) = A^[l + {qv*)-P]ve-^f\ (D4) 

with A = 0.322, p = 0.3, and q = 0.707. The density threshold v is set to S c /a(M, z) with S c = 1.686. Then, the halo 
biasing parameters can be calculated from Eq. (|D4j) as 

b\{M,z) = l + ex+Ex, (D5) 



where we define 



%(M,z) = — (e 1 + E 1 )+e 2 + E 2 , (D6) 



qv 2 - 1 qv 2 qv 2 - 3 
(\ = — 7 , (2 = — x j , (D7) 

Or. Or Or 



2p 1 E 2 _ 1 + 2p 

Xl + {qv 2 )P' ~E 1 = ~ 

In the expression (|D3[I . there appears the minimum halo mass, M m j n , which can be determined from the condition, 



dM n h {M,z){N) M . (D9) 

Af mi „ 

Finally, it seems rather difficult to determine the fiducial value of the remaining parameter N , because physical 
meaning of the parameter TV is less clear. In this paper, we just adopt the relation (|C5|) . and compute N assuming 
^0 = 0: 



N =\j w? p ^ 2 - (D10) 

APPENDIX E: SYSTEMATIC BIAS FOR THE BEST-FIT PARAMETERS 

In this appendix, we briefly review how to estimate the biases in best-fit parameters arising from the systematic 
effects. We arc especially concerned with the impact of neglecting massive neutrinos on the dark energy constraints. 
In this case, the biased parameter estimation is obtained by fitting the observational data to the power spectrum 
template incorrectly assuming f v = 0. Let us write down the observed power spectrum as 

P° hs (k) = P^°(k) + P™ isc (k) 

= Pl"=°(k)+PI* s (k) + P™ ise (k), (El) 

where P^ olse (k) denotes the instrumental noise, while the systematics in power spectrum, P| ys (fc) is defined as 
P| ys (fc) = P^"^°(fc) — Pf" =0 (k). Then, the systematic bias in a certain parameter Sp a is computed in (e.g., (95[) as 

S Pa = ^(F- 1 )^, (E2) 





d 3 q 
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galaxy 
a/3 



F™ B ) in which the neutrino parameter /„ is excluded 



where F is the full Fisher matrix (namely, F Q ^ = F 

from the matrix element. Note that in computing F, the fiducial parameter for f v must be set to /„ = 0, because 
we consider the situation that the observed power spectrum is incorrectly fitted to the template neglecting massive 
neutrinos. Here, the vector quantity S a is represented S% MB + Sg alaxy , which are respectively given by 



S CMB = ^Y^Cf^i^C^Cj)}-^, 

i x,y Pa 
V s {z t ) 



ggalaxy _ 



47T 2 



( Zi ) pcst,sy S(fc . ^ 91n pcst (fc . J 

k z dk 



dp a 



n g (zi)P?*(k;zi 



(E3) 
(E4) 



where the angular power spectrum for CMB, Cf' sys , is defined similarly to the case of galaxy power spectrum, i.e., 
c x, sys _ c x,f„ji0 _ C x,f„=Q^ Note again that we get f v = Q i n computing Cf and P„ est . 



[1] E. Komatsu et al . larXiv:0803. 05471 [astro-ph], 

[2] S. Dodelson, Modern Cosmology, Academic Press, San Diego (2003). 

[3] S. Fukudaet al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 85, 3999 (2000). 

[4] S. N. Ahmed et al. [SNO Collaboration], Phys. Rev. Lett. 92, 181301 (2004). 

[5] K. Eguchi et al. [KamLAND Collaboration], Phys. Rev. Lett. 90, 021802 (2003). 

[6] T. Arakki et al. [KamLAND Collaboration], Phys. Rev. Lett. 94 081801 (2005). 

[7] R. D. Mc Keown and P. Vogel, Phy s. Rept. 394, 315 (2004). 

[8] B. Kayser, |arXiv:hep-ph/0506165| 

[9] J. Bonn et al., Nucl.Phys.B (Proc. Suppl.) 91, 273 (2001). 
[10] K. Ichikawa, M. Fukugita and M. Kawasaki, Phys. Rev. D 71, 043001 (2005). 
[11] M. Fukugita, K. Ichikawa, M. Kawasaki and O. Lahav, Phys. Rev. D 74, 027302 (2006). 
[12] M. Kowalski et al, Astrophys. J. 686, 749 (2008). 

[13] D. J. Eisenstein et al, Astrophys. J. 633, 560 (2005); W. J. Percival et al, Astrophys. J. 65 7, 51 (2007); W. J. Percival 
et al, Mon. Not. Ro y. Astron. Soc. 3 81, 1053 (2007); E. Gaztanaga, A. Cabre and L. Hui. larXiv:0807.3551l [astro-ph]; 
W. J. Percival et al. larXiv:0907.1660l [astro-ph.CO]. 

[14] J. R. Bond, G. Efstathiou and J. Silk, Phys. Rev. Lett 45, 1980 (1980). 

[15] W. Hu, D. J. Eisenstein and M. Tegmark, Phys. Rev. Lett. 80, 5255 (1998). 

[16] 0. Elgar0y et al, Phys. Rev. Lett 89, 061301 (2002). 

[17] M. Tegmark et al, Phys. Rev. D 74, 123507 (2006). 

[18] K. Ichiki, M. Takada and T. Takahashi, Phys. Rev. D 79, 023520 (2009). 

[19] T. D. Kitching et al, Phys. Rev. D 77, 103008 (2008). 

[20] U. Seljak, A. Slosar and P. McDonald, JCAP 0610, 014 (2006). 

[21] J. R. Pritchard and E. Pierpaoli, Phy s. Rev. D 78, 065009 (2008). 

[22] M. Sumiyoshi et al, larXiv:0902. 20641 [astro-ph.CO], 

[23] M. Takada, E. Komatsu and T. Futamase, Phys. Rev. D 73, 083520 (2006). 
[24] S. Hannestad and Y. Y. Y. Wong, JCAP 0707, 004 (2007). 
[25] F. B. Abdalla and S. Rawlings, Mon. Not. R. Astron. Soc. 381, 1313 (2007). 
[26] D. Jeong and E. Komatsu, Astrophys. J. 651, 619 (2006). 

[27] R. Angulo, C. M. Baugh, C. S. Frenk and C. G. Lacey, Mon. Not. Roy. Astron. Soc. 383, 755 (2008). 

[28] R. Takahashi et al, Mon. Not. Roy. Astron. Soc. 389, 1675 (2 008). 

[29] A. G. Sanchez, C. M. Baugh and R. Angulo. [arXiv:0804.023"3l [astro-ph], 

[30] H. J. Seo, E. R. Siegel, D. J. Eisenstein and M. White. larXiv:0805. 01171 [astro-ph], 

[31] D. Jeong and E. Komatsu, arXiv:0805.2632 [astro-ph]. 

[32] T. Nishimichi et al, arXiv:0810. 0813] [astro-ph]. 

[33] P. McDonald, Phys. Rev. D 75, 043514 (2007). 

[34] P. Valageas, Astron. Astrophys. 465, 725 (2007). 

[35] S. Matarrese and M. Pietroni, JCAP 0706, 026 (2007). 

[36] T. Nishimichi et al, Publ. Astr on. Soc. Jap. 59, 1049 (2007). 

[37] M. Crocce and R. Scoccimarro, larXiv:0704.2783l [astro-ph] . 

[38] T. Matsubara, Phys. Rev. D 77, 063530 (2008). 

[39] A. Taruya a nd T. Hiramatsu, Astrophys. J. 674, 617 (2008). 

[40] M. Pietroni. larXiv:0806.097T1 [astro-ph], 

[41] H. Nomura, K. Yamamoto and T. Nishimich i, JCAP 0810 031 (2008). 

[42] J. Carlson, M. White and N. Padmanabhan. larXiv:09 05.0479 [ astro-ph .CO]. 

[43] A. Taruya, T. Nishimichi, S. Saito and T. Hiramatsu. larXiv:090 6.0507l [astro-ph.CO], 



32 



J. Lesgourgues and S. Pastor, Phys. Rept. 429, 307 (2006). 

F. Bernardeau, S. Colombi, E. Gaztanaga and R. Scoccimarro, Phys. Rept. 367, 1 (2002). 

S. Saito, M. Takada and A. Taruya, Phys. Rev. Lett. 100, 191301 (2008). 

P. McDonald, Phys. Rev. D 74, 103512 (2006). 

M. Tegmark et al, Astrophys. J. 606, 702 (2004). 

A. Taruya, Astrophys. J. , 537, 37 (2000). 

S. Pueblas and R. Scoccimarro, arXiv:0809.4606 [astro-ph]. 

C. P. Ma and E. Bertschinger, Astrophys. J. 455, 7 (1995). 
U. Seljak and M. Zaldarriaga, ApJ, 1996, 469, 437 

A. Lewis, A. Challinor and A. Lasenby, Astrophys. J. 538, 473 (2000). 

D. J. Eisenstein and W. Hu, Astrophys. J. 511, 5 (1999). 
W. Hu and D. J. Eisenstein, Astrophys. J. 498, 497 (1998). 
F. Bernardea u, Astrophys. J. 4 33, 1 (1994). 

R. Takahashi. [arXiv:0806. 14371 [astro-ph] . 

T. Hiramatsu and A. Taruya, Phys. Rev. D 79, 103526 (2009). 
N. Makino, M. Sasaki, and Y. Suto, Phys. Rev. D 46, 585 (1992). 

B. Jain and E. Bertschinger, Astrophys. J. 431, 495 (1994). 
Y. Y. Y. Wong, JCAP 0810, 035 (2008). 

A. Kiakotou, O. Elgaroy and O. Lahav, Phys. Rev. D 77, 063005 (2008). 

R. E. Smith et al. [The Virgo Consortium Collaboration], Mon. Not. Roy. Astron. Soc. 341, 1311 (2003). 

J. Brandbyge, S. Hannestad, T. Haugboelle and B . Thomsen, JCAP 0808, 020 (2008). 

J. Brandbyge and S. Hannestad, arXiv:0812.3149 [astro-ph]. 

R. Takahashi et al., Astrophys. J. 700, 479 (2009). 

J. N. Fry and E. Gaztanaga, Astrophys. J. 425, 1 (1994). 

J. N. Fry, Astrophys. J. 461, L65 (1996). 

M. Tegmark and P. J. E. Peebles, Astrophys. J. 500, L79 (1998) 
L. Hui and K. P. Parfrey, Phys. Rev. D 77, 043527 (2008). 

A. F. Heavens, S. Matarrese and L. Verde, Mon. Not. Roy. Astron. Soc. 301, 797 (1998). 
U. Seljak, Mon. Not. Roy. Astron. Soc. 318, 203 (2000) 

R. E. Smith, R. Scoccimarro and R. K. Sheth, Phys. Rev. D 75, 063512 (2007) 

C. Alcock and B. Paczynski, Nature, 281, 358 (1979). 

H. J. Seo and D. J. Eisenstein, Astrophys. J. 598, 720 (2003). 

T. Matsubara and A. S. Szalay, Phys. Rev. Lett. 90, 021302 (2003). 

L. Perotto, J. Lesgourgues, S. Hannestad, H. Tu and Y. Y. Y. Wong, JCAP 0610, 013 (2006). 
M. Tegmark, Phys. Rev. Lett. 79, 3806 (1997). 

A. Klypin, J. Holtzman, J. Primack and E. Regos, Astrophys. J. 416, 1 (1993). 
M. Takada and B. Jain. larXiv:0810.4170l [astro-ph], 

E. Sefusatti, M. Crocce , S. Pueblas and R. Scoc cimarro, Phys. Rev. D 74, 3522 (2007). 
[Planck Collaboration], arXiv:astro-ph/0604069 

T. Matsubara, Astrophys. J. , 615, 573 (2004). 

N. Kaiser, Astrophys . J. 284, L9 (1984). 

A. J. Albrecht et al, |arXiv:astro-ph/0609591| 

J. Lesgourgues, S. Matarrese, M. Pietroni and A. Riotto, JCAP 0906, 017 (2009). 
A. Lewis and A. Challinor, Phys. Rept. 429, 1 (2006). 

D. Munshi, P. Valageas, L. Van Waerbeke and A. Heavens, Phys. Rept. 462, 67 (2008). 
S. Singh and C. P. Ma, Phys. Rev. D 67, 023506 (2003). 

A. Ringwald and Y. Y. Y. Wong, JCAP 0412, 005 (2004). 
K. Abazajian et al, Phys. Rev. D 71, 04 3507 (2005). 

M. Shoji and E. Komatsu. larXiv:0903.2669l [astro-ph.CO] . 

R. Scoccimarro, R. K. Sheth, L. Hui and B. Jain, Astrophys. J. 546, 20 (2001). 
R. K. Sheth and G. Tormen, Mon. Not. Roy. A stron. Soc. 308, 119 (1999). 

B. Joachimi and P. Schneider, arXiv:0905.0393 [ astro-ph.CO]. 



|http://www.science. doe.gov/hep/hepap/feb2007/hepap_bennett-feb07.pdf 
http: / /wigglez. swin.edu.au/ 
http: / /cosmology.lbl.gov/BOSS / 

http://arxiv.org/ftp/astro-ph/papers/0510/0510272.pdf 
http:/ /www.as. utexas.edu/hetdex/ 

http: / / sci.esa.int /science-e /www / object / index.cfm?fobjectid=42266 
http:/ /www.science. doe.gov/hep/hepap/feb2007/ hepap_bennett_feb07.pdf 

To compute the results shown in Fig. [T] we used the fitting formula of the linear growth function, D c h(k; t), developed in 
Eq.(12) of [55|, for computational simplicity. 
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FIG. 11: The 2nd-order and 3rd-order growth functions. The ratios of 2nd-order growth functions to the square of 
linear-oder growth functions, Ag(ki, k2)/(D c b(ki)D c b(k2) (top left) and Bg(k\, k2)/ (D c b(ki)D c b(k2) (top right) are plot- 
ted (blue curved surfaces). In the case of the 3rd-order growth functions, the ratios to the cubed linear-oder growth functions, 



Tg / (D c b(ki) D c b(k2)), (T = A — F) are shown (blue curved surfaces). Although the 3rd-order growth functions, T s depend 
on three specific wavenumbers, it is sufficient to specify two wavenumbers to calculate P^ 3 (k). That is why we show the 

ratios such as the combination of k\ and k2, Tg 3 \ki, &2, fci) and lf\k\, k\, kz). Note that T s 3 \k\, k2,k\) = lf\k\, fci, fo) for 
X = B, D and F. For the reference, the constant values (red plane for 2nd-order and green plane for 3rd-order) are also shown. 
These constant values corresponds to the SPT treatment in a CDM model. 



